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ACOUSTICAL  SCATTERING  FROM  THE  IMPEDANCE  COVERED  STRAIGHT  EDGE 
AND  WEDGE;  THE  EXACT  THEORY 

ABSTRACT 

The  analytic  extension  of  the  Fourier  integral  from  the  half  space 
to  vedge-like  spaces  leads  to  Sommer feld  type  integrals.  Theorems  derived 
by  Malyuzhinets  in  conjunction  with  new  methods  of  handling  Fourier  Laplace 
integrals  make  it  possible  to  apply  the  boundary  conditions  to  the  integrand 
to  determine  the  integrand  function  and  thus  to  determine  the  solution 
from  the  boundary  conditions.  The  resulting  Malyuzhinets  functions  can 
then  be  closely  approximated  by  simple  expressions  for  real  and  complex 
arguments.  The  diffracted  field  then  depends  on  the  shadow  boundary  and, 
in  contrast  to  the  classical  theory,  also  on  the  orientation  of  the 
diffracting  surface.  The  diffraction  field  is  described  by  the  sum  of 
two  terms.  One  term  can  be  interpreted  as  a  scattering  effect  because 
of  the  discontinuity  in  the  medium  caused  by  the  impedance  surfaces;  the 
second  could  be  interpreted  as  a  field  generated  by  reemission  of  some 
of  the  energy  that  is  bent  into  the  diffractor  near  the  edge.  This  term 
is  equal  to  half  the  reflected  amplitude  at  the  shadow  boundary  of  the 
reflected,  wave  and  to  half  the  amplitude  of  the  transmitted  wave  at  the 
shadow -boundary  of  the  transmitted  wave .  It  describes  the  trivial  condition 
of  continuity  in  the  transition  fiom  the  insonr-ified  region  to  the  shadowed 
region.  On  the  side  of  the  incident  wave,  the  amplitude  is  very  nearly 
equal  to  the  excitation  as  though  due  to  the  reflection  factor  that  would 
obtain  if  the  angle  of  exit  (of  diffraction)  were  produced  by  reflecting 
a  ray  on  the  we dge  surface  under  the  appropriate  angle  of  incidence.  As 
a  consequence,  this  second  term  vanishes  »t  two  angles  (the  Brewster 
angles)  that  are  equidistant  from  the  normal  of  the  surface.  In  this 


report,  the  diffraction  is  investigated  of  a  great  number  of'  straight 
edges  and  wedges  whose  surfaces  are  covered  with  various  impedances. 
The  exact  solutions  are  compared  with  approximations.  However,  even 
the  exact  solution  can  be  easily  derived  because  of  a  simple  and  very 
good  approximation  to  the  Malyuzhinets  functions. 


I.  INTRODUCTION 


The  Fourier  Integral  for  the  Half  Space  and  Its  Analytic  Extension  to  Wedge 


Spaces 

In  dealing  with  propagation  problems,  we  solved  the  wave  equation  in  a 
suitable  coordinate,  system.  Because  the  wave  equation  is  of  second  order, 
we  obtain  two  independent  solutions.  Boundary  conditions  are  usually  prescribed 


,4.  -  - - 


coordinate  surfaces  and  at  infinity*  One  solution  frequently  lias-  to  te 


uiscarded  because  it  ispresents  energy  sources  at  infinity.  In  the  case  of  a 
straight  edge  (semi  infinite  plane ),  boundary  conditions  need  to  be  prescribed 
at  both  surfaces,  i.e,  the  solution  must  be  discontinuous  at  the  half  plane. 
The  coordinate  surface  is  no  longer  a  plane  hut  is  a  degenerate  cylindrical 
paraboloid. 

The  standard  Fourier  method  is  not  applicable  to  wedge  spaces  nor  to 
the  straight  edge.  Fourier  integrals  comerge  only  in  the  semi  space. 

The  Fourier  integral  can  be  built  up  from  plane  two-dimensional  waves: 


p  =  S  e 
o 


j(wx+  k^y)  -  gut 


(1.1) 


where  S  is  an  amplitude  constant,  k  is  the  wave  number,  and 
o 


->  /  ,  2  2 
K  =  V  XL  -  W 

y 


The  quantity  w  can  be  interpreted  as  the  forced  wave  number  in  the  x 
direction  at  the  surface  y  =  0;  the  quantity  k  is  the  wave  number  of  the 
wave  that  propagates  in  che  medium  as  a  consequence  of  the  vibration  forced 
with  the  wave  number  w  along  the  plane  y  -  0.  Wo  may  now  prescribe  the 
vibration  amplitude  at  y  =  0  by  a  Fourier  integral  and  extend  the  integral 


to  the  y-space  by  including  the  k  part  in  the  exponent.  Tims  we  have 
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p  „  /•  sj(«<  +  yVk!i  .  w2)  3o(w)du  _  (i.3) 

-00 

where  Sq(w)  is  the  spectral  amplitude  function.  The  integral  then  represents 
a  solution  of  the  wave  equation  provided  that  it  is  absolutely  convergent. 

This  is  the  case  in  the  upper  half  plane  if  the  exponent  has  a  negative  real 
part  for  y  >  0.  and  if  S  (w)  is  bounded  in  the  interval  of  integration. 

The  convergence  properties  of  the  integrand  are  not  affected  by  a 
coordinate  transformation.  However,  it  will  be  of  advantage  to  introduce 
cylindrical  coordinates  (see  Fig.  1.1): 

x  =  r  sin  cp  ,  y  =  r  cos  cp  ,  w  =  k  sin  a  ,  (l.l) 

Regardless  of  the  sign  of  the  square  root  in  the  integrand  of  Eq.  (3), 

the  transformation  introduces  once  more  an  ambiginity  of  the  sign 
(Vr  -  sin~a  =  +  cos  a).  The  exponent  in.  the  integrand  thus  becomes 

jkr  (sin  a  sin  cp  +  cos  cp  cos  a)  =  +  jkr  cos  (cp  +  a)  (1.5) 

The  limits  are  given  by 

w/k  =  sin  a  =  sin  (a  +  ja.  )  =  sin  a  cosh  a.  ,i  cos  a  sinh  a.  =  +•  «  , 

(1.6) 

We  find  that  cos  a  =  0  ,  a.  =  H  k>  ,  a  =  $  for  w  =  «  ,  and  that 
r  ’  1  —  r  2  ’ 

n  ,  ,, 

a  =  -  —  ,  a.  =  +  a,  for  w  =  -  00. 

r  2  x  — 


We  have  a  choice  of  two  basically  different  values  for  each  limit. 


It  does  not  make  much  difference  which  sign  we  select  for  the  square  root; 
convergence  can  be  obtained  for  each  sign  by  selecting  properly  the  new 
limits  of  integration.  If  we  select  the  positive  sign  for  the  square  root, 
the  integral  will  converge  if 


kr  sin  (cp  -  a  )  sinh  a. 

v  r  i 

(1.7) 

is  positive.  This  means  that  we  must  have 

0  <  cp  -  a  <jt,  a.  — =■,  w 
r  x 

(1.8) 

-  n  <  cp  -  a  <0,  a.— >  -  o> 

r  i 

(1.9) 

For  cn — 4  oo  we  obtain  for  the  limit-n/2  of  cp: 

0  <  -  n/2  -  <  n  or  -  ~  <  ar  <  -  n/2  , 

(1.10) 

and  for  the  limit  n/2  of  cp: 

0  <  n/2  -a  <n  or  -  ^  <  a  <n/2 

i*  o  r 

(1.11) 

Thus  a  =  -  ~  -  c  when  cp  is  negative  (c— >u)  and  u  n  -  —  +  when  cp  is 
r  u  t  Cj 

positive,  or  we  may  assume  that  ar  =  ~  ^  ?  an<^  that  9  is  restricted  to  the 
interval  -  (--  -  e)  <  cp  <  (i  -  c ) .  Similarly  we  find  that  -  +  n/2  for 
a.  — >  -  oo.  The  transformed  integral  is  thus  given  by 

11  j  - 

-  2  + 

P(r,cp)  =  e3l<:i  cos^  a ^  SQ(k  sin  a)  k  cos  a  da  (1.12) 


irt 
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If  we  assume  that  the  path  of  integration  passes  from  +  .■  ~  to  -  — 

and.  from  there  along  to  -  n/2  and  then  straight  up  to  -  ^  then  the 

contribution  of  the  path  along  the  real  axis  represents  the  field  of  waves 
in  the  direction  cp  that  propagates  at  an  angle  -  ~  <  a  <  ~  with  respect 
to  the  real  axis.  The  contribution  of  the  vertical  parts  of  the  path 
represents  the  field  of  waves  in  the  direction  cp  that  propagates  along 
the  axis  a.  ~  0  and  decays  strongly  with  distance.  We  have  no  freedom  in  the 
path  of  integration  if  the  integral  is  to  converge  for  all  values  of  cp  that 
correspond  to  point.-,  in  the  upper  half  plane.  Thus  changing  the  variable 
did  not  change  anything  basically. 

Equations  (1.1)  and  (1.5)  derive  the  field  from  the  boundary  condition 
at  y  -  0.  This  boundary  condition  is  equivalent  to  a  field  that  is  generated 
by  a  source  distribution  somewhere  between  y  =  0  and  y  =  -  <».  In  a  wedge 
space  or  in  case  of  a  straight  edge,  we  need  solutions  that  converge  also 
for  y  <  0  and  that  are  discontinuous  at  the  surfaces  of  the  straight  edge 
or  wedge.  This  requirement  excludes  solutions  in  Cartesian  coordinates. 

Such  solutions  have  the  period  2 it  and  would  lead  to  the  same  values  at  the 
two  surfaces  of  the  straight  edge.  If  we  want  to  perform  an  analysis  in 
terms  of  simple  harmonic  functions,  .,hey  must  have  a  period  that  differs 
from  2 it ,  V/e  are  therefore  certain  that  B  (a)  4  S  (a  +  2it). 

To  introduce  the  preceding  condition  explicitly  into  the  solution, 
we  write  the  integrand  in  the  following  form : 

S(a  +  it)  -  S(a  -  m)  =  2itj  SQ(k  sin  a)  k  cos  a  ,  (1.13 

S(a  -  n)  can  then  he  interpreted  as  the  radiation  from  the  front,  S(a+  n) 
from  the  back  (see  Fig.  1,2).  Front  and  back  are  separated  by  the  serni 
infinite  plane,  each  point  therefore  contributes  because  it  radiates  into 


the  front  space  and  because  it  also  radiates  into  the  shadow  space.  This 
particular  Geometry  seems  to  be  impressed  on  the  whole  space.  The  physical 
space  is  defined  ha.?;  as  the  space  -  n  <  q>  <  n  (Sonmierfeld) .  All  other 
angles  (;T|  >  it )  represent  a  purely  mathematical  space  (the  Riemann  space 
in  the  Sommerfeld  theory,..  The  reflected  fields  can  then  be  interpreted  as 
being  generated  eg  mirror  image  sources  that  are  hidden  in  the  Riernan  space ; 

this  agrees  with  the  assumption  of  having  only  one  source  in  the  physical 
space . 

The  solution  then  is  represented  by  the  expression 


where  we  have  introduced  the  new  variable  a'  =  a  +  it  ,  a"  =  a  -  x  in  the  two 
parts  of  the  integral,  written  a  again  for  a'  and  q",  and  changed  the 
limits  of  integration  correspondingly.  The  new  paths  of  integration  are 
shown  in  Fig.  1.3a.  We  now  assume  that  5(a)  is  regular  in  the  range  of 
integration  except  for  a  number  of  discrete  poles  at  finite  distance  from 
the  real  axis.  We  combine  the  two  halves  of  the  path  in  each  half  plane 
into  a  loop,  getting  the  path  ?'(<p)  as  shown  in  Fig.  1.3b.  Poles  at  infinite 
distances  from  the  real  axis  can  be  acc.o\mtea  for  by  considering  their 
residues.  We  thus  obtain: 
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It  is  inconvenient  to  introduce  limits  that  depend  on  cp.  By  changing 
variables 


(cp  —  a)  —  -  u  j 


u  =  a  -  cp 


a  =  u  h-  cp  ,  da  =  du 


(1.16) 


and  writing  a  again  for  u  the  integral  transforms  into 


/  -jkr  cos  a 
/  e 

r 


5  (a  +  cp)  du 


(1.17) 


we  are  free  to  move  the  path  of  integration  within  the  regions  of  convergence 
of  the  integrand.  Since  the  convergence  of  the  integrand  does  not  depend  on  cp  , 
the  path  ?  is  the  same  for  all  values  of  -  0  <  cp  <  b  .  We  have  assumed 
implicity  that  S(a  +  cp)  vanishes  sufficiently  strong  at  infinity  so  that  we 
are  allowed  to  start  the  path  at  +  and  to  proceed  in  any  direction  as 
long  as  we  stay  in  the  convergence  range  of  the  integrand.  For  positive  , 
these  limits  are  now  given  by 


kr  sin  a.  sinh  a.  >0,  0  <  a  <  ji  , 

l  l  ’  r  5 


(1.18) 


if  a.  >0  and  by 

-  2x  <  a  <-jt,  if  a.  <0 

r  5  i 

An  analogous  result  is  obtained  f  i  the  negative  half  plane.  The  path  I' 
is  as  shown  in  figure  3.. 4.  It  extends  along  the  real  axis  from  at  least 
-  it  to  +  it.  It  is  apparent  that  the  function  S(a)  must  be  analytic,  for 
at  least  the  range  -  ix  <  (a  +  G>)  <  «  for  ac  ">  c  =  const., 


7 


or 


-  (jt  +  C1)  <  ar  <  it  +  0 


(1.19) 


(i.e.  except  for  the  poles  near  the  real  axis.  We  thus  have  arrived  at  a 
generalized  Fourier  integral  which  will  converge  in  wedge-like  spaces.  We 
shall  refer  to  this  integral  as  the  loop  integral  in  the  following. 


8 


II.  PROPERTIES  OF  THE  LOOP  INTEGRAL 
(a)  Derivative  with  Respect  to  ep 

The  derivation  of  the  loop  integral  is  best  performed  in  its  form  (1.14 ) . 
Since  a  small  displacement  of  the  path  y(cp)  into  7(q>  +  drp)  within  the  con¬ 
vergence  range  of  the  integral  does  not  change  its  value,  differentiation 
can  be  confined  to  that  of  the  integrand.  Thus 


ol 

QCP 


P  5  e“Jkr  COS  ^ 

J  Sp 

7(9) 


s(a)da 


(2.1) 


Z'  -jkr  cos 


jkr  e  ^  —  -  sin  a  +  cp)da 


(2.2) 


where  we  have  replaced  cp  -  a  by  -  a’  and  written  a  for  a'  again. 


(b )  Condition  for  Zero  Value  of  the  Integral 

The  most  important  theorem  that  applies  to  the  loop  integrals  is 
Malyu.zhi net's  theorem  I.  A  loop  integral  is  zero 

/-jlr  cos  z  /  \  ,  „  .  „ . 

e  0  s(z)dz  =  C  (2.5) 

r 

if  the  function  s(z)  is  even,  i.e.  if 

s(z)  =  o(-z)  (2,4  ) 

where 


z  =  u  +  jv 


To  prove  this  theorem,  the  path  of  integration  is  deformed  into  the  two 
lines  infinitely  close  to  v  =  -it  and  u  =  -Hr  (see  Fig.  l.Jj-)  and  integration  is 


performed  either  to  the  left  or  to  the  right  of  these  lines  so  that  the 
exponent  has  a  negative  real  part.  Since  cos  (z)  =  cos  (-z)  and  s(z)  -  s(-z), 
the  integrands  along  the  line  0  >  v  >  -op,  u  =  -  x  and  u  =  it,  .0  <  v  <  »  are 
equal.  But  since  the  directions  of  the  two  parts  of  the  path  are  opposite, 
their  contributions  cancel. 

If  the  function  s(z)  increases  at  a  high  rate  towards  infinity,  condition 

■X- 

(5)  has  to  be  replaced  by  a  more  complex  condition. 


( c )  Inversion  Formula,  for  the  Wedge  Integral  for  cp  =  const . 

Malyuzhinets  also  derived  an  inversion  formula  for  the  wedge  integral. 
If 


.  2  /  y.  ~  it  \ 

is  -  ni  s  =  0  -  <  arg  in  <  — ) 


(2.5) 


_ ,  y  1  /  m r  cos  a  „ ,  x , 

S(r)  =  2^-r  J  e  S(a)da 


(2.6) 


then 


=  - 


m  sin  a  D  „ ,  \  -nr  cos  a.  . 

- 2 - J  DVr;  U 


/  r\  rr  \ 

•  1 ) 


The  derivation  is  straightforward  and  given  in  the  reference:  G.  D.  Malyuzhinets 
Inversion  Formula  for  the  Sommer f eld  Integral,  Mathematical  Soviet  Physics 
Melody  3;  1958,  92-5*1 . 
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III.  BOUNDARY  CONDITIONS 

The  function  S(a)  is  determined  from  the  boundary  conditions  for  the 
sound  pressure  at  the  surfaces  cp  =  +  <t> .  For  a  wedge  with  pressure  release 
surfaces 


P(r ,4)  =  y>rcos  a  S(a  +  9)  da  =  0  (3.1) 

r 

To  satisfy  this  relation,  S(a  +  $)  must  be  an  even  function  of  a.  Thus 

S(a  +  <t>)  s=  S(-a  +  $) 

S(a  -  4>)  =  S(-a  -  <t>)  (3.2) 

If  we  replace  a  by  a  -  $  in  the  first  and  by  a  +  4>  in  the  second  equation, 
the  left  hand  sides  become  equal  and 


S(a)  =  S(-a  -  2<b)  -  S(-a  +  2®) 


(3.3) 


This  result  shows  that  the  solution  has  the  period  4$.  As  pointed  out  in 
Section  I,  Equation  1.19,  the  solution  must  be  analytic  in  the  band 
|Re(a)|  <  'V  +  jr  except  for  a  pole  for  cp  =  cp^,  whose-  contribution  represents 
the  incident  wave.  The  simplest  function  that  has  a  period  k't>,  has  a  pole 
at  a  =  cp  .  Such  a  function  is 


(a)  = 


sin  —  (a 


(3.4) 
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But  this  function  does  not  satisfy  condition  (3.3)  because 

S(-a  +  ?.<!>)  =  sin[“  (-a  +  20  -  <pQ ) ]  =  sin[^  (-a  -  cpQ)+  «]  =  sin  ~  (a  +  cpQ)  ^  S(a) 

(3.5) 

But  the  product  of  two  functions  of  the  type  S(«),  one  with  +  cp^,  the  other 
with  -  cp^  in  the  argument  will  satisfy  our  requirement  because  of  the  symmetry 
which  is  thus  generated.  We  are  thus  led  to  assume  that 


S 


(a)  = 


sin  ~  (a  -  <pt 


,)  sin^-  (a  +  <Po) 


(3.6 


\ 

) 


S(a)  still  has  a  pole  at  a  =  -  cp^;  but  q>  =  -  (pQ  represents  an  incident 
wave  in  the  physical  space  (-  $  <  <p  <  $).  But  there  is  no  second  incident 
wave,  and  this  pole  must  be  cancelled  by  multiplying  by  a  periodic 
non-trivia!  factor,  A,  that  has  such  a  pole.  The  factor 


A  =  sln  W,  «  +  sin  ~  90 


(3.7) 


fulfills  all  requirements.  It  cancels  the  pole  and  satisfies  the  conditions. 
Hence  we  have 


S(u)  ~ 


B(sin  ~  a  +  sit,  ^  cpQ) 

Gltl  4  (a  +  Gin  ib  (a  -  <(,o) 


(3.8) 


The  residue  of  L>(a)  at  a  -  qj  is  unity  if 


r,  n  Jt 

B  =  z5T  C0G  -  % 


(3.9) 
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We  thus  arrive  at  the  solution 


s(a)  =  ~ 


cos  cpc 


It 
il  0 


sin  —  (a  -  cp  )  cos  L_  (a  +  cpQ) 
li  0  °  7T5 


cos  .  (a  -  cp  )  cos  (a  +  cp  )  +  sin  *  (a  -  cp  )  sin  ^  (a  +  cp  ) 
),Q _ °  )[Q _ ° _ H  °  10 _ 

sin  —2-  (a  -  cp  )  cos  LL  (a  +  cp  ) 
to  0  to  ° 


(3.10) 


-li-  [  cot  JL.  (a  -  cp  )  +  tan  _7C  ( a  +  cp  )] 

H  1$  °  40 


For  the  stationary  phase  integration,  S(a)  is  "best  written  in  the  form 


S(a)  = 


cos  cp 
21'  'o 


•pth 


yJ  .  J.,1.  j 


sin  ~  a  ~  sin  ^  cpQ 


which  is  easily  deduced  from  the  third  form  on  the  right.  The  cosine  factor 
in  the  numerator  is  a  consequence  of  the  pressure  release  boundary.  This 
factor  is  zero  for  cp^  =  0  .  We  thus  have  obtained  the  well  known  Sommer  fold 
integrand  for  the  pressure  release  wedge;  in  construct  to  Sommer fell’s 
heuristi  procedure  the  method  above  is  based  only  on  mathematical  deduction. 
The  solution  for  the  rigid  wedge  can  be  derived  in  a  very  similar  manner. 


:  -1 


n 

cos  n 
20 

*o 

[sin  n  a  +  sin  n  cp 
20  20  ° 

] 

,,  i 

i 

W 

sin  _L_ 
20 

(a 

+  qj)  sin  LL  (a  -  cp  ) 

°  20  ° 

r 

n 

cos  (— 2— 
20 

V 

1  2  sin  (a  +  cp  )  cos 

4  0  0  ’ 

JL  (a  - 
5o 

CPo) 

i 

2<J> 

t  sin  — — 
to 

-  (a  +  cp  )  cos  n  (a  4-  cp  ) 

0  j,*  0 

sin  Tf 
to 

(a  -  cp  )  cos 

! 

LL  (a  -  cp  )  ! 

il  o  0  1 

•  esk.u 
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Malyuzhinets  has  derived,  a  new  kind  of  Fourier-Laplace  analysis  which 
can  also  be  used  to  derive  the  function  S(a).  But  in  case  of  the  Sommer- 
feld  wedge,  the  Malyuzhinets  Method  is  impractical  because  of  the  complexity 
of  the  integrals  that  result  and  the  preceding  method  is  preferable. 


1/1 


IV.  THE  FOURIER  LAPLACE  ANALYSIS  OF  MALYUZHINETS 

Functions  that  increase  exponentially  towards  +  oo  or  approach  finite 
values  at  +  co  can  be  represented  by  Laplace  integrals  by  separating  them 
into  two  parts,  one  part  representing  the  function  in  the  interval  0  to  co, 
the  other  zero  part  representing  the  function  in  the  interval  -  »  to  0. 

This  procedure  leads  to  two  different  spectral  functions .  To  compute 
inverse  transform  because  of  their  different  regions  of  convergence,  each 
of  these  two  functions  has  to  be  integrated  over  a  different  path.  The 
analysis  that  has  been  introduced  by  Malyuzhinets  is  based  on  a  single 
spectral  function  F(w)  for  the  whole  range  of  integration.  However,  the 
representation  of  this  function  depends  on  whether  the  line  integral  is 
at  the  left  or  right  of  the  imaginary  axis.  This  simplification  becomes 
possible  because  the  fields  we  are  interested  in  either  approach  finite 
values  o i  vanish  at  ^  uu.  The  integrals  that  occur  in  the  various  computa¬ 
tions  usually  contain  sines,  cosines,  tangents,  and  groups  of  functions 
of  a  variable  z  that  are  regular  in  strips,  a  <  Re  (z)  <  b.  For  instance, 
l/cos  z  is  regular  for  — jt/2  <  Re  (z)  <  ^  .  We  therefore  limit  our  study 
to  functions  which  exponentially  approach  constant  values  for  Im(z)-^  -i-« 
when  a  <  Re  z  <  b.  We  thus  have  by  definition 

f(z)— f(jcn)  H-  o  {e  ( -  )  j  „  f(j00)  +  0<e^Z',  Im(z )  -x  (4.1) 

+  0  (e-^Iin(z)  }  =  f(-jr»)  +  0  (e"duz)  ,  lm(z)~  -fa  (4.2) 


X,  p  positive  constants 
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The  analysis  we  are  going  to  develop  could  be  extended  to  exponentially 
increasing  functions  as  j  Xm ( z )  |  — >  co.  However,  we  would  then  loose  the 
effect  of  the  constants  f(j»)  and  f(-jco).  That  wave  fields  differ  at  + 
infinity  is  easily  illustrated  for  a  half  plane  as  diffractor,  half  the 
incident  energy  will  be  reflected  back  at.  the  semi  infinite  plane.  But  at 
great  distances  from  the  plane,  the  shadow  space  will  be  filled  up  with 
energy  by  diffraction  and  the  amplitude  will  approach  l/h/2  times  that  of 
the  incident  wave.  Thus  the  wave  amplitude  will  be  at  +  »,  and  S 
at  -  « , 

The  function  f(z)  -  f(-j»)  can  be  represented  by  a  two-sided  Laplace 
integral  because  it  vanishes  exponentially  for  Im(z)-^-  «,  and  approaches 
a  constant  for  z  -  +  300.  To  derive  this  integral,  let  us  use  the  notation 
of  Malyuzhinets  and  prove  that  the  integral 


dz 


F  (w)  =  P  [f(zj  -  f(-jco)]  el 

converges  for  0  <  Re(w)  <  p  ,  and  that  the  integral 


0  <  Re  w  <  u 


(4.3) 


F  (w)  =  ~~p=r  p[f{z)  -  f(jw)]  e^WZdz  ,  ~7\  <  Re(w)  <  0 

Vj>  tt  is. 

—  Too 


(4.4) 


converges  for  -1  <  Re(w)  <  0.  To  perform  this  proof,  we  divide  the  integral 
(4.3)  into  the  two  integrals 


j 


s 

-.1“ 


[f(z)  -  f(-j»)]eJWZdz  +  — 

ViT 


J" 

f 


[f(z)  -  f(-joo)]  e^WZdz 


(4 


The  first  integral  exists  because  f(z)  is  analytic  and  f(z)  -  f(-o°°)  vanishes 
at  its  lower  limit.  If  we  replace  f(z)  by  Eq.  (4.1)  the  first  integral 
reduces  to 


o 

p  [f(-»  -  f(-»  +  0  {e"jUZ]  ]  ejv2dz 
O' 

-0“ 


This  integral  converges,  provided  R e(w)  =  p  <  0.  In  the  second  integral  we 
write  correspondingly 

Pi  f(jco)  -  f(-»  +0  (e^z)]  GJwzdz  (4 

U 

o 

This  integral  converges  if  Re  w  >  0  and  Re(w)  +  A  >  0.  The  three  conditions 
then  lead  to 


0  <  Re  w  <  p 


The  preceding  formulae  can  be  inverted  as  is  proved  by  methods  well 
Xnown  in  Fourier  analysis 

„**/  V  -.iwz  . 

•c  v  w )  e  ~  aw  ,  (  m 

provided  -  A  <  Re(w)  <  0  and 


■Pf  rr 


■Pf  A..  \ 


■i  n 

•vtt  J. 


-.1° 


f(z) 


f(-» 


0°° 

— P  F*(w)  e~^WZ  dw 
-D°° 


provided  0  <  P.e(w)  <  |i 


If  we  subs  tract  the  two  equations  from  each  other- 


17 


d®+e  j„-e 

ty1  F*(w)e-J«d*  .  />  F"(„)e-J“V) 

-joo4^  ^  • 


(1.11) 


-jm-e 


where  we  have  replaced  the  condition  Re(w)  >  -  >,  by  shifting  the  path  of 
integration  by  e  to  the  left  and  correspondingly  the  condition  Re(w)  <  n 
by  shifting  it  by  e  to  the  right  where  e  >  0  car.  „■  a-s  small  as  we  please. 
Since  f(z)  is  analytic  in  the  strip  a  <  Re(z)  <  b,  F*(w)  and  F**  (w)  cannot 
have  any  singularities  other  than  for  Re(w)  =  0,  and  cannot  have  singularities 
for  Im(w)  f  o  since  the  integrand  is  analytic  in  the  strip  and  finite  at 
infinity;  the  only  admissible  singularity  is  a  simple  pole  at  w  =  0  since 
F(w)  is  single  valued.  We  have  for  w  /  0,  by  applying  Eq.  (4.3), 


F(w)  <  litn 


dz  =  lim 


finite 


(4.12) 


The  right  hand  side  of  Eq.  (4.11)  reduces  to  the  integral  around  the  pole 
at  w  =  0  and  F(w)  becomes  infinite  for  w  =  0.  The  residue  of  F(w)  at  this 
pole  is  given  by 


residue  • 


f (-j“)l 


(;..13) 


^  -x-  -x-* 

tor  all  values  oi  z.  mis  snows  snai-  r  —  v  .A  more  rigorous  proof  can 
be  derived  with  the  aid  of  the  Cauchy  Chilbert  transform  representation  of 
the  Fourier  transform. 

Vie  can  still  improve  the  symmetry  of  the  reverse  transform  by  adding 
Eqs.  (4.9)  and  (4.10)  and  dividing  by  two: 


f(z)  =  i  [  f(j«)  +  f(-D»)]  +  — 4=  [  P  F  (w)e  u  dw 

2  'dY2«  U  . 

-joo-e 


jco+C 

-jm+e 


F  (w)e  *^WZdw] 


(4.14) 


where  F*(w)  and  F**(w)  are  given  by  Eqs.  (4.3)  and  (4.4).  We  can  drop 
the  star  and  the  double  star.  F(w)  is  an  analytic  function  which  in  the 
band  -  A  <  Re(w)  <  0  can  be  represented  by  the  integral  F**(w)  and  in  the 

-X- 

band  0  <  Re(w)  <  p  by  the  integral  F  (w).  Equation  (4 .14)  represents  a 
very  convenient  and  also  the  most  general  form  of  the  solution.  But  for 
the  final  evaluation,  the  two  integrals  have  to  be  combined  and  written  in 
a  standard  form.  This  can  be  done,  for  instance,  by  shifting  the  path  of 
integration  of  the  second  integral  to  the  left,  and  adding  the  contribution 
of  the  residue  at  w  =  0  to  the  result  (see  Fig.  4.1): 


-e+jco 

f(z)  =  r  Cf(ooo)  +-  f(-j®)]  +  —j= 

2  2  V2n 


+  •  residue  at  (w  =  0)  , 

%  Ct  Jl 


[F*(w )  +  F**(w)] 


-1WZ 


dw 


(4.15) 


or  by  changing  w  into  -  w  in  the  second  integral  of  Eq.  (4.14): 


-e+joo 

f(z>  +f(-a->]  +4=  ,P  [F*(»)e'i“adw  4-  F«(-W)^wba»  (1.16) 


Frequently,  integration  can  be  further  simplified  by  discarding  the  odd  parts 
of  the  integrals  but  retaining  possible  contributions  at  the  pole  w  =  0. 
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V.  THE  IMPEDANCE  WEDGE 


(a)  The  Solution 

In  this  section,  wedges  will  be  investigated  whose  surfaces  are 
covered  by  acoustic  impedances.  The  reflection  factor  for  a  plane  incident 
wave  at  a  plane  surface  of  acoustic  impedance  z,  is  give^  by 


—  COS  CD  +  1 

pc 


n  Pc 
cos  9  - 

2 

PC 

cos  cp  +  — 


cos  cp  -  cos  0 
COS  cp  -r  COS  0 


(5.1) 


where  cp  is  the  angle  of  incidence  with  the  normal  to  the  surface.  For 
the  angle  given  by 


COS 


cos  -  -S)  =  sin  i3  , 


(5.2) 


the  reflection  is  zero.  The  angle  d  =  it/2  -  0  is  defined  as  the  grazing 
acoustic  Brewster  angle,  0  as  the  Brewster  angle  with  the  normal  to  the 
surface.  This  definition  applied  formally  also  for  complex  impedances. 

In  deriving  the  acoustic  equations  for  wedges,  it  is  of  considerable  advantage 
to  measure  angles  from  the  central  wedge  plane  (see  Fig.  1.1),  The  wedge 
surfaces  then  are  given  by  cp  =  +  0  .  The  boundary  conditions 
that  have  to  be  satisfied  are:  (I)  For  the  upper  surface,  where  cp  = 


v 

9 


_P_  _ 1_  . 

z+  ~  kpc  d(rcp) 


(time  factor  e  ; 


or  if  we  introduce  the  Brewster  angle 


(5-3) 


iwaMfrpa mtmmvu  maamm 
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jk  p  sin  d+  -  i  (gj)^ 


For  the  ..  ower  surface  where  cp  =  -  <P,  z  =  z_;  -6  =  t3_ 


jkp  sin^  =  -  i  (^)<p=_* 


(5.4) 


(5.5) 


If  we  apply  these  boundary  conditions  to  the  integral  solution  and  using 
Eqs.  (2.1)  and  (2.3),  the  following  equations  result: 

(sin  a  +  sin  t3+)  S  (a  +  0)  =  -  (sin  a  -  sin  -3+)  S  (  -a  +  4>)  =  0 

(sin  a  +  sin  -3_ )  S  (a  -  0)  =  -  (sin  a  -  sin  -S_ )  S  (-a  -  G>)  =  0  (5-6) 


It  is  apparent  from  the  Kirchoff  theory  that  the  impedance  cover  of  the 
faces  of  the  wedge  will,  not  displace  the  shadow  boundary  of  the  wedge  surfaces 
significantly  and  it  is  obvious  that  it  cannot  effect  the  pole  nor  its  residue 
that  represents  the  incident  wave.  If,  therefore,  we  write  the  solution  in 
for  form 


s(a)  =  £l£)l(c£) 

¥(<P0) 


(5-7) 


where 


a(cr)  is  the  Sommerfeld  form  of  the  integrand  of  the  loop  integral 
for  a  pressure  release  wedge, 

']/ (a)  is  a  function  that  describes  the  effect  of  the  impedance  cover, 
i|/(cp  )  is  a  kind  of  normalizing  factor, 

then  o(a)  will  contain  all  the  poorly  behaving  parts  of  the  solution  (poles 
describing  incident  and  reflected  waves)  and  ~  -y  will  be  a  well  behaving 
function  with  no  poles  within  trie  range  of  integration  (but  it  may  have  poles 
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at  the  extremes  a  =  +  it).  We  shall  see  in  section  (  XII)  that  the  solution 
of  the  impedance  wedge  must  be  based  on  the  pressure  release  wedge  and  not 
on  the  rigid  wedge.  The  impedance  wedge  backscatters  at  all  angles.  The 
rigid  straight  edge  does  not  backscatter  when  the  wave  impinges  parallel 
to  the  plane  of  the  straight  edge  so  that  cp  =  0.  The  resultant  solution 
would,  therefore,  break  down  for  cp^  =  0  if  cr(o:)  were  the  integrand  for  the 
rigid  wedge.  If  we  enter  Eq. (5-7) into  the  two  boundary  equations  above, 
we  obtain  two  equations  of  the  form  of  Eq.(5-6)for  the  function  i|r(a),  S(a) 
being  replaced  by  i(a).  The  cr(a)  part  cancels  out  left  and  right  because 
0\a)  satisfies  the  conditions  (3.2)  for  a  pressure  release  wedge.  To  solve 
equation  (5 .6),  we  take  the  logarithm  left  and  right  and  differentiate. 
Because  the  logarithm  of  (-1)  is  jrt,  and  its  derivative  is  zero,  we  obtain: 


cos  a  :  +  _4>) _ cos  « _ V  (-a  +  ^ 

sin  a  +  sin  \p(a  4-  )  sin  a  -  sin  '  \|r(-a  + 


(5-8) 


A  similar  equation  results  with  -3+  replaced  by  i8_  and  <P  replaced  by  -  <1*. 
Next  we  introduce  the  new  function: 


(5-9) 


The  Malyuzhi nets  transform  of  f(a  +  'I')  and  f(-tt  +  <I>)  then  are  given  by 


( Eq .  h  .3  anu  ,h  ) : 


Mf(a  +  <1>)  =  e  ^W'  F(w) 


Mf(-a  -i  i<)  =  e'"1  J  F(-w) 


(5-10) 
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where  mf(a)  -  F(w)  is  the  Malyuzhinets  transform  of  f(a)  and  the  letter  M 
means  symbolically  Malyuzhinets  transform  of  f(a).  Let 


„  (a)  _ _ £os  a  _  cos  a _ 

sin  a  +  sin  -8  sin  a  -  sin  -8 


2  sin  9  cos  a  _  k  sin  9  cos  a 
sin2 a  -  sin2 -9  =  cos  2a  "  Co~ 


2  sin  9  cos  a 


-  sin(a  +  -9)  sin  (a  -  9)  sin  (a  +  9)  sin  (a  -  9) 


(5.11) 


where  9  is  the  grazing  Brewster  angle  and  0  the  Brewster  angle  with 
respect  to  the  normal  to  the  surface.  We  have  omitted  the  bars  above  the 
symbols  9  and  6  that  mark  them  as  complex  quantities.  The  transformed 
equations  (5-8)  then  are 

g+(w)  +  F(w )  +  F(-w)  =  0 

g_(w)  +  F(w)  +  F(-w)  --=  0  (5.12) 

where  g+(w),  g_(w)  represent  g(w )  -  Mg  (a)  with  9  replaced  by  9,,  9 
respectively.  Elimination  of  F(-w)  leads  to: 


6+(w)  e 
F(w)  =  - 


■jwl  ,  jwd- 
-  g  (w)  e 


2j  sin  2w'l 


(5.13) 


Because  e(-0")  and  g(+j'0  are  zero,  the  inverse  transform  of  the  first  term 
of  Eq .  (5.13)  then  is  represented  by 


f(z)  = 


3'" 

-h 

fi  it  J 
-3"' 


da 


n 


•> 

^■3”’  j 


L.' 


-£-Jr" 


■:~3'- 


__ _ 1  .  c;iW^'Z"'^dv 

in(^”H-  9)  ~  sin  (a  -  9)  2j  sin(~2f!v)- 


S  H  3 


(5.  HO 
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where  we  have  written  for  i9+.  In  the  second  integral  we  replace  w  by 
-w.  The  two  paths  then  coincide,  and  the  two  contour  integrals  can  be  combined: 


f(z)  = 


J00  0°° 

P  t  p  da  [7^  1 


^  J  u 
-j00  -  j°° 


'sxnla  -i- 


“57 


_ 1  i  sin  w(a-z-$)  .  . 

sin  (a  -  d)  sin  (2i*w)  dw  (5*15) 


Next  we  replace  a  +  3  by  a  +  d  -  |  =  a'  [a  =  a'  +  (|  -  ,  in  the  first  factor, 

and  a  -  d  by  a  -  -S  +  ~  =  a"  ,  [a  =  a”  -  (t-  -  ■3)]  in  the  second  factor. 


f  ( z  )  =  - 


1 

57 


/  y7 

-j®  -j“> 


da’ 


cos  a 


r  sin  w[a'  -  (z  4-  $  +  d  -  ■-)  ]  — 


dw 


2  sin  2Cw 


J" 

7/  / 3in  w[a"  -(*  +  *-  *  +  |)1 
-d~  -O'" 


dw 


<s  sm  2iw 


0-16) 


Since  there  is  no  pole  at  -jt/2  >  a  <  jc/2  ,  we  were  allowed  to  discard  the 
real  parts  of  the  limits  of  integration.  We  can  perform  the  v  integration 
first.  The  following  integral  is  of  a  tabulated  form: 


>- 

f')  sin  va  . 

/  — - —  aw 

■j  mi  wb 


sin  hvb 
o 


dv^j  r  tan 


2b 


(5.17) 


Performing  the  v  integration  in  the  first  integral  Ey.,(p.l6)  we  obtain: 


cos  a 


„  ,  n  1  ft  <  yO.  -(z-t-'M'1  -  n/2)n  1 

*  -  Sr  J  tan[ - si - 1  — d 

-0" 


where 


i3  ‘bfz 


tc  n  v  57 

cog"  "( v~|7; 


-  dv  (7 • ) 


-jc 


\  1 


2k 


a-(z  +  <l,  +  i3-^-)=-v 

a  =  -  v  +  ji 

n=z+a+d-|  (5.19) 


The  w  integration  in  the  second  integral  leads  to  a  function  f  (z)  which 
differs  only  in  that  S  -n/Z  is  replaced  by  ~ ( 1  -  %/Z).  Two  more  such 
functions  then  would  result  for  ($— i  -$),  the  face  <p  =  -  which  are 
similar  in  form  except  that  3  — >  d- . 

Four  integrals  are  of  different  form,  they  are  obtained  if  we  perform 
the  a  integration  first.  In  this  case  it  is  expedient  to  combine  the  two 
integrals,  Eq.  (5.1 6)  and  to  write  Z  =  z  -1-  C>: 


f(z) 


4* 


99 


do*  ^sin  v* f g:  — 


a 

“  > 


i°L 


„  *7 


cos  a 


sm 


iW 


+  c'' 


dw 


-Joo  -30 


1  P P sin  w(a  -  Z)  cos  w£  dw 

2jc  t/</  cos  a  sin  2lw 

-.1“ 


jco  jco 

1  P  /)( sin  mo.  cos  w  Z  -  cos  wa  sin  wZ)cos  wO  dw  do: 
fit  ij tj  cos  a  sin  2tw 

*0“  -j00 


cos  wa  sin  wZ  cos  wO 
cos  a  sin  2^w 


dadw 


(5.20) 


because  only  the  even  part  of  the  a  integral  contributes  to  the  integral. 
The  a  integration  is  of  a  tabulated  form: 


1 J 


n 

/ 


wa 


da 


cosh  wu 
cosh  u 


du 


cos  a 


cos 
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and 


in  w Z  cos  w0 


cos  w  —  sin  2$w 
J“  2 


dw 


-ft 


sinh  sZ  cosh  s$ 


ds 


cosh  s  —  sinh  2$s 
o  2 


(5.21) 


_  _  1  p _ sinh  s(z  +  0) 


cosh  (s  $)  sinh  24>s 
o  c 


1  f)  sinh  s(z  -  0)s 

OS  “  rr  / - > - - -  dS 

2  ^  cosh  (s  |)  sinh  2is 


(5-22) 


=  f(z{)  +  f(z') 


where 


«  o 

z|  =  Z  +  0  =  z  +  <I>  +  %  '  ® 

z'  =  Z  -  0  =  z  +  4>  -  (£  -  -9)  (5-23) 

2  <£ 


and 


f(z)  =  -  i 


sinh  zsds 


(5-2*0 


^  ^  cosh  s  §  sin  24>e 
o  2 


Again,  two  more  terms  result  for  the  contributions  of  the  face  <p  -  -  O 
which  are  similar  in  form  except  that  4>  is  replaced  by  -  *,  and  t  denotes 
d_.  The  logarithmic  derivative  of  the  solution  thus  consists  of  the  sum  of 

four  terms 


d  lot;  »(z)  =  t'f-l  =  f(z 
dz  *(z) 


+  H  0  )if(z,  -1-  4-  -  0+)+f(z  -  1  + 


0  )-l-f(z  -<!>  -  0_) 

(5-25) 


■n-Tnir - ■— . m  ■  mmimtm *&*&**■**»*, 

\  ‘  *  f 
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We  still  have  to  integrate  each  term  up  to  the  desired  value  of  z  .  We 
assume  arbitrarily  z'  =  0  for  the  lower  limit  of  integration.  This  pro¬ 
cedure  is  admissible  because  the  resultant  solution  is  represented  by 
4'(a)/vy (cp  )  so  that  the  integration  constants,  which  lead  to  identical  factors 
in  v(a)  and  i|/(cP0)  cancel  out  in  the  final  expressions.  We  thus  have 


where 


and 


(5.26) 


(5.27) 


because  of  the  exponential  nature  of  the  solution;  the  exponents  add,  the 
functions  multiply.  We  may  still  write  the  arguments  in  full: 


'i-(a)  =  i',Jso.  +  Q  +  |  -  fl+)  ^(a  +  <b  -  |  +  -8  )  -  4>  +  |  -  d  )  ^(a  -  *  -  |  i  *  ) 


+  *  +  0+)  ^(a  +  -  o+)  t4(a  -  4  +  Qj  \if  (a  -  4-  -  o_) 


(5-28) 


The  solution  of  the  diffraction  problem  then  is  given  by 


?(r,T,k) 


1_  P  -ikr  cos  (a 
2ni  U  G 
r 


q'^S(a)da 


(7.29) 


where 


'■Hot) 

Vlvc ) 

U 


S(a)  =  a(«) 
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it© 

'  o 

COS  — 7— 

20 


_  f  .  it  ,  s  x  «(«  h-  cp  )]  ^  (5.30) 

^  (a  -  <PQ)  +  tan— ^ 


£in  if  ■ sin  zi 


The  exact  solution  is  given  by 


P(r,©)  =  l 

r 


1  /J-jkr  cos  a  v/(a  +  cp)  rr(a  +  cp) 


,  U  ^  j  o  /  .  nu  -r 

tt?T  oos  (m - p- 


TO?  _1 

sin  -— -)  da  ,  (5.31) 


If,  assuming  kr  >  0  ,  we  deform  the  integration  contour  1  into  two  paths  of 
steepest  descent  through  the  saddle  points  a  =  +  it. 


P(r/P)  =  |J  £ 


«  (■*  .  f\.  \ 

j  i^r  ~t-  51/ h  / 


V  2ickr 


i  'H©  -  it)  .  „  t(©  +  it)  ! 

i'T  '  nip-  z  ~~<^r  j 


(5.32; 


where  and  M  are  the  Malyuzhinets  angle  factors 


it (q)  -  it) 

sm  ^  _  s:L1 


(5-33) 


.  it  (tp  H-  it )  .  11  *o 

sm  -1-1 — - -  -  sin  t- 

20  20 


(5.34) 
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VI.  THE  MALYUZIUHETS  \|u(a)  FUNCTIONS  AND  THEIR  APPROXIMATION 

The  Malyuzhinets  function  ty(ce)  is  the  solution,  Eq.  (5 .6)  without  poles 
or  zeroes  in  the  strip  Re!  a!  <  4>.  They  are  usually  computed  by  developing 
the  integral  solution  into  a  series  up  to  the  pole  of  the  integrand,  or  by 
evaluating  the  integral  directly  on  the  computer.  Beyond  the  pole  of  the 
integrand,  the  integral  no  longer  represents  the  solution  of  the  problem, 
and  one  of  the  extension  formulae  (see  A.ppendix  A)  must  be  used  to  determine 
the  functions  for  greater  values  of  the  agreement.  Figure  (6.1)  shows  ^(z) 
for  real  z  =  a  +  <3>  +  0+  as  a  function  of  z;  Fig.  (6.2)  represents  the  same 
function  ^  (z)  and  Fig.  (6.3)  the  phase  angle  for  complex  argument  as 

a  function  of  the  real  part  with  the  imaginary  part  x  of  z[Jm(z)  as 
parameter ] . 

The  tyj,  functions  are  relatively  insensitive  to  the  argument  a  except 
for  argument  values  near  their  zeroes  and  their  poles.  A  small  imaginary 
part  in  the  surface  admittance  increases  the  magnitude  at  the  zero  to  a 
finite  small  value,  and  reduces  the  magnitude  at  the  pole  to  a  finite  value; 
the  effect  on  the  resulting  function  is  practically  negligible  if  Im  pc/z  < 
(see  Fig.  6.2).  'Die  surface  then  acts  as  if  its  impedance  were  real.  If 
the  imaginary  part  of  the  relative  surface  admittance  pc/z"  exceeds  1,  the 
curves  become  almost  horizontal  and  the  wedge  acts  as  if  it  were  pressure 
release.  Because  of  this  relatively  small  sensitivity  of  the  diffraction 
field  on  the  impedance  on  the  surface,  we  can  expect  that  the  results  also 
apply  to  boundaries  that  cannot  be  described  by  point  impedances,  such  as 
plate  boundaries. 

Computations  with  the  Malyuzhi nets  functions  are  greatly  simplified 
by  the  relations  that  exist  between  them,  and  which  are  summarized  in 


rol  j—1 


appendix  A.  Furthermore ,  there  is  no  need  to  evaluate  the  defining 
integrals.  Since  the  Malyuzhinets  functions  are  smooth  functions,  they 
can  be  approximated  to  a  high  degree  of  accuracy  by  an  expression  of 
the  type 


,v. 

COS 


a 


2  a 


■) 


COS 


a 


Va> = 


zero 


/* 

cos  (r- 


Gt 


(6. 


2  a 


pole 


where  a  and  a  ^  are  given  by  Eq,  (A22).  This  function  reproduces  the 
zero  pole 


exact  slope  and  magnitude  of  i^(a)  for  a  =  0,  and  is  zero  at  a  and  infinite 


at  the  pole  at  the  second  cosine  factor  in  the  numerator  then  is  used 


to  improve  the  fit  in  the  pole  region.  We  find  the  following  values  for 


0.3* 


l/3rt 


0.4* 


3/8* 


0.5* 


0.7n  0.9* 


4.63 


5.04 


5.87 


5.56  6.98 


8.6 


9.72 


10A9 


These  functions  approximate  i|r  (a)  in  the  whole  range,  also  beyond  the  pole 
and  in  the  complex  domain.  Fig.  (6.4)  and  (6.5)  show  the  functions  $n(°) 
and  4'  ,  (0)  as  computed  by  Eq.  (6.1).  The  exact  values  have  been 

31  /  C 

computed  by  Zavadskii  and  Sakharova  which  fall  on  the  two  curves. 
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VII.  THE  ANGLE  FACTORS  OF  THE  MALYUZHINETS  SOLUTION 


The  angle  factor  in  the  Malyuzhinets  integrand  results  by  adding  the 
two  angle  factors  that  occur  in  the  Sommerfeld  integrand  for  the  pressure 
release  wedge  (see  Eq.  3.10).  The  Malyuzhinets  stationary  phase  solution 
then  contains  the  angle  factors 


-  cos  — r  cp 

^  =  - - - 

sin  §r  (cp  -  *)  -  sin  ~  «PQ 


cos  rr- 

M,  - - - 

sm  —  (cp  +  n)  -  sin  ^  9q 


(7.1) 


We  shall  assume  cp^  is  positive.  For  negative  cp^,  the  two  angle  factors 
simply  interchange.  The  angle  factor  is  plotted  in  Fig.  (7.1)  for  a 
straight  edge  and  in  Fig.  (7*2)  for  various  wedges  with  the  angle  of 
incident  as  parameter.  If  4>  >  a/2,  has  no  poles  in  physical  space 
-  T  <  cp  <  <3>  because  the  first  term  in  the  denominator  then  is  always 
negative.  increases  towards  the  limits  cp  =  +  4>  of  the  physical  space 
without  reaching  a  maximum  (where  the  slope  is  horizontal).  The  two 
greatest  values  of  then  are  given  by 


“  cos  21 


■1-4)  n 

-  +  cos  n 


3in  h 


(7-2) 


then  has  a  minimum  for  (cp  -  n)n/ 2$  =  -jt/2  or  cp  =  n  ~  when  the 
diffracted  ray  propagates  along  the  continuation  of  the  face  cp  =  -<I>  of 


the  wedge, 


^1  min 


cos  f? 

1  +  sin  h  % 


<  003  20  % 


(7-3) 


and  then  increases  constantly  towards  the  faces  cp  -  +  0  of  the  wedge  to  the 
value 


"  COS  20  *0 

±  COS  "  sin  cpo 


(7» 


The  second  angle  factor 


_ COS  2~  *0 

(cp  +  if)  -  sin  ~  <p 


(?.?) 


has  two  polos,  one  at  the  shadow  boundary  of  the  reflected  wave  (cp  -  20 

•*  n  -  rpo)  and  one  at  the  shadow  boundary  of  the  incident  wave  (cp  -  -  jc  -1-  cp  ). 

It  has  a  true  minimum  between  the  two  shadow  boundaries  for  cp.  $  _  h: 


cos  2*  % 

(M J  .  =  - - - — 

x  2  'man  .  .it 

1  -  sun  2Z  <P0 


and  it  decreases  from  the  shadow  boundaries  towards  the  wedge  at  cp  =  -1-  <f 
to  the  same  values  as  those  given  by  Eq.  (7. 2)  for  at  +  4>. 
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If  cpQ  is  negative  and  cp  is  replaced  by  -  cp,  changes  into  and 
into  M^.  Fig,  (7.1a)  shows  the  basic  shape  of  the  curves  M^(rp)  and 
Mg('!')  and  the  angular  ranges  in  which  each  term  is  positive  and  negative. 

Figure  (7, Id)  to  g  show  for  comparison  also  the  Sommerfeld  and  the 
Kirchhoff  solutions  for  the  rigid  and  the  pressure  release  straight  edge. 

We  find  that  is  not  very  different  from  the  Kirchhoff  or  Sommerfeld 
solution  for  the  pressure  release  straight  edge.  Figures  (7-3)  and  (7 A) 
represent  the  angle  factors  for  backscatter  when  cp  =  qiQ  as  rectangular  and 
polar  plots,  M0  and  turn  out  to  be  the  exact  mirror  images  of  each 
other  at  the  axis  cp  =  0.  We  can  confine  our  attention  to  the  space  0  <  cp  < 
because  when  looking  at  the  backscattered  field  we  can  always  assume  we 
are  in  front  of  the  plate  (<p  >0).  The  maximum  of  the  angle  function 

M0  then  represents  the  diffraction  at  the  shadow  boundary  of  the  reflected 
wave  when  it  is  coincident  with  the  incident  wave,  j .e.  when  it  impinges 
normal  to  a  wedge  surface.  This  is  the  case  when  the  transducer  is  just 
at  the  shadow  boundary  of  the  geometrically  reflected  wave.  In  practical 
situations,  the  transducer  will  be  within  the  cone  of  the  reflected  waves 
or  it  will  be  outside.  Inside,  diffraction  does  not  contribute  much  to 
backscatter.  Outside  of  it,  the  received  signal  is  exclusively  due  to 
diffraction.  The  amplitude  of  the  received  signal  then  decreases  greatly 
with  the  angular  distance  from  the  shadow  boundary. 

For  pressure  release  surfaces: 

pc/z  =  cos  0  =  cos  (0  +  jA  )  =  cos  0r  cosh  0_.  -  j  sin  0 ^  sin  G^_-a,ou  (7,7) 

Because  the  real  part  of  the  impedance  must  be  positive,  0_^  <  —  ,  we  must  have 

0 .  —  *4*  00 

1  — 


But  as  <x-Aw,  'J/^(ar+  jo.  )  — 7 .  ( jot  )— ■ (see  Fig,  6.2  )  and  the  functions 

in  the  numerator  and  denominator  cancel.  Thus  the  multiplier  T(a  + 


35 


and  the  solution  is  -,iven  by 


_rt  e  .i(kr  +  £) 

2'v  it  2rtkr” 


+  Mg) 


(7-8) 


The  two  angle  factors  then  simply  add  up,  and  the  resultant  expression  reduces 
to  the  classical  solution  for  a  pressure  release  wedge. 

For  a  rigid  wedge  =  0_  =  0  ,  and  by  Eq.  All 


*(a)  =  *t(f)U  [cos  &  (a  +  *)  cos  (:*  -  *)Z 


1  _  jt  It  * 

-  2  Lcos  a  +  cos  ^  * 


:•]  - 


1  /  It  \  ^ 

2  Vg*  cea  H  L< 


(7-9) 


ind 


^  (a)  cos 


■  3<h 


v7c?  ;  "  n 
v  To  cos 


(7.10) 


For  a  straight  edge,  the  functions  f  (cp  -  it)  and  \|;(cp  +  it)  then  become 
proportional  to  sin  cp/2  and  to  -  sin  cp/2,  and  the  two  terms  in  the 
Malyuzhinets  solution  counteract  each  other  in  the  angular  range  cp  >  0  . 
Thus,  the  solution  for  the  rigid  straight  edge  becomes 


P  = 


-jt  e 


2'1‘ 


■(Jkr  +  g) 
V  2rtkr 


sin  cp/2  sin  cp/2  Mg 

cos  cp  Tz  Cos  cp  /2  \ 
o'  o' 


(7.11) 


Comparison  with  the  corresponding  Sornmerfeld  solution  gives  the  Sommer f eld 
factors  in  terms  of  M-j.  and  Mg,  (see  Appendix  B).  The  converse  procedure, 
i.e.  the  expressing  of  the  factor  and  Mg  in  terms  of  the  Sornmerfeld 


**MW‘ 


\ 
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factors  leads  to  very  complex  expressions  that  contain  also  the  angle 

p 

functions  of  (rt/2<l')(a)  aud  of  n L‘jz$. 

For  an  impedance  matched  wedge : 

sin  $  =  1  ,  -3  =  n/2  ,  6  =  0  (7.12) 

+  ■’  i  + 

and  by  Eq.  A13 


f(a)  =  4^,(0)^  ^y2(a)  (7'13) 

The  t  function  then  is  proportional  to  the  square  of  the  function 

that  is  plotted  in  Fig.  5.1a.  For  a  straight  edge,  20/2  =  «  and  the  zero 
of  this  function  occurs  for  a  =(3/2)n  ,  and  that  of  the  solution  for 
cp  +  -n  =(3/2 )rt,  9  =  it/2;  i.e,,  there  is  no  field  reflected  or  diffracted 


normal  to  the  surface  of  the  straight  edge. 
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VIII,  THE  MALYUZHIHETS  FUNCTION  AND  THE  CLASSICAL  REFLECTION 
FACTOR  AND  THE  EXIT  FACTOR 

The  second  term  In  the  Malyuzhinets  stationary  phase  solution  (Eq,  5-32) 
is  closely  related  to  the  classical  reflexion  factor.  If  cpQ  is  the  angle  of 
incidence  with  respect  to  the  plane  of  symmetry  of  the  wedge,  the  angle 

of  the  incident  ware  with  the  normal  of  the  wedge  surface(cp  =  <I>  -  ji/2) 
is  given  by  (see  Fig.  8.1): 

<PNo  "  *  "  */2  -  <PQ  (8.1) 

and  the  angle  of  the  reflected  wave  with  the  cp  =  0  axis  by 

cpR  =  4>  -  n/2  +  cpNo  =  2i>  -  «  ~  <po  (8.2) 

For  this  angle,  the  angle  factor  M  (Eq.5.34),  in  the  Malyuzhinets  solution 
Eq,  (5.32)  has  a  pole  which  determines  the  field  at  the  shadow  boundary. 

To  evaluate  the  solution  for  the  shadow  boundary,  the  integral  (5.31 )  is 
written  in  the  form 


V 


n(s) 

Z  ~  Z 


rg(z), 
e  dz 


(8.3) 


The  stationary  phase  solution  then,  when  z  coincides  with  the  stationary 
phase  point,  is  given  by 


P 


jit  h 


,  v  rg(z  ) 
( z  )  e  J  ■  o 


(8.M 


(see  Jones,  The  Theory  of  Electromagnetism,  Pergamon  Press  iy6h ,  New  York, 
p.  690).  Thus  we  find  from  (5.31)  for  the  contribution  of  the  saddle  point 
Re(z)  =  +  it: 


3  6 


,  jkr  ( it  +  9 )  cos(jt<p Jz<b)  j 

P(r,tp)  =  Jn  - 


¥T 


jt  it  /  ,  > 

24  COS  24  (*  4  Cp) 


(8.5) 


<P- 


R 


=  lA 


,3kr  -  %,» 

+(<v 


* 

It  is  shovm  in  the  classical  theory  of  diffraction  that  this  term  should  be 
equal  to  half  the  amplitude  of  and  of  opposite  signs  as  that  of  the  reflected 
wave . 

Thus,  we  are  justified  to  interpret  the  ratio 


-  y(24  -  T0) 


R(mo) 


(8.6) 


as  the  reflection  factor  of  the  wedge  surface  for  a  wave  incident  under  an 

-■  * 

angle  (4  -  cp^)  with  respect  to  the  plane  of  the  reflecting  surface.  To 
prove  this  conclusion,  let  us  assume  for  reasons  of  simplicity  that  both 
wedge  surfaces  are  covered  with  the  same  material  so  that  Q^=  0__  =  0  . 

We  then  have  (see  Eqs .  5.28  and  A13 ) : 

♦(c)  h/7g  j  Tj)  V«<B  - 

AAA  1I1  .  (m  +  O')  ■ii .  ,„(r I;  -  0 ) 

,4/2',0  '  W/c,  o 


(8.7) 


and  using  Equation  AI5 : 


Akudrzyk :  Foundations  of  Acoustics,  page  570-571*  Springer  Verlag,  Dev/  York 
and  Vienna,  1972. 


tmrnm 


to.  a  '**■*"*  vA" 
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¥(2$  -  9q) 
V(fPQ) 


y2(24)  -  +  0)  y2(2*  -  - 0) 

\/o(-^  +  77 )  %./■><■ 


<172 


=  cot  1/2  (-cp  +  0  +  It/2  +  $)  cot  1/2 (-cpo  -  6  +  jt/2  +  ®) 

cos  0  -  sin  (-  cp  +  $) 

_ _ _ 2 - -  =  .  B(rp  ) 

cos  0  +  sin  (-  cp  +  ‘M  c 


(8.8) 


But  i  -  cpQ  is  the  grazing  angle  of  the  incident  wave  with  the  reflecting 
surface  of  the  wedge.  If  we  replaced.  i>  -  cpQ  by  the  angle  with  respect  to 
the  normal  of  the  surface,  the  sine  would  be  replaced  by  the  cosine  of 
chat  angle.  Thus,  R(cpQ)  is  identical  with  the  classical  1  eflection  factor 
of  the  impedance  covered  infinite  surface.  This  classical  reflection  factor, 
therefore,  describes  the  effect  of  diffraction  at  the  shadow  boundary. 

The  impedance  covers  of  the  wedge  surfaces  modify  the  diffraction 
phenomenon.  This  modification,  relative  to  the  diffraction  of  a  pressure 
release  wedge,  is  described,  by  the  functions 


ty(cp  +  it) 


and 


(cp  -  it) 

V(t|>o) 


(8.9) 


Because  y(cp  )  is  a  function  that  varies  very  little  with  the  angle  of  incidence, 
o 

* 

the  above  two  functions  are  practically  independent  of  the  angle  of  incidence. 
This  conclusion  is  the  most  remarkable  result  derived  from  the  exact  theory. 

It  is  the  angle  cp  of  exit  of  the  diffracted  rays  that  determines  the  in¬ 
fluence  of  the  impedance,  regardless  of  the  angle  of  the  incident  radiation 
(see  Fig.  8.2). 


The  factor  cos  jfj>  /2l  occurs  also  in  the  classical  theory  if  the  field  is 
formulated  in  the  Maly uzhinets  form. 
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For  an  arbitrary  exit  angle  cp,  we  express  the  angle  of  exit  by  what 
would  be  the  angle  cp^  of  incidence  for  cp  as  an  angle  of  exit  in  the 
second  term  in  the  solution;  i.e.  cp  =  24  -  ji  -  cp^.  The  second  term  of 
the  solution  then  takes  the  form 


ty(cp  H-  it) 


*(2*  -  cp’)  ^(cp^) 

V(<p;r  F(qT7 


(8.10) 


It  is  proportional  to  the  factor  R(cp^)  which  we  shall  define  as  the  exit 
factor.  The  factor  in  the  rectangular  bracket  is  very  nearly  unity  as 
long  as  the  angle  of  incidence  cp  and  the  angle  of  incidence  cp'  (if  the 
diffracted  ray  were  a  reflected  ray)  are  somewhat  different  from  +  4. 

Figures  8.3  show  a  comparison  of  the  classical  exit  factor  R(cp^)  with 
the  Maiyuzhinets  function  tJ/ (cp  +  n).  If  we  divided  -Jc (cp  +  jr)  by  *(m  )  it 
becomes  exactly  equal  to  R(cp^).  If  wc  move  the  Maiyuzhinets  curve  up  or 
dovm  in  the  logarithmic  plot,  the  two  curves  coincide  for  cp  =  cp^,  and  then 
the  difference  between  the  two  curves  will  give  us  the  error  that  results 
if  we  replace  ^ (cp  +  n ) /'[< (qM  )  by  R(q/).  The  curves  prove  that  this  error 
is  practically  negligible  everywhere  except  in  the  region  near  cp  =  0.  We 
thus  obtain  a  very  good  approximation  to  the  second  term  by  replacing 

M ,'k (cp  +  x)/'KrP,.)— *M,,R(cp* ) 

in  the  range  q>  >0,  cp^  >  0.  The  impedance  cover  reduces  the  diffracted 
ray  by  the  exact  factor  which  is  the  same  as  the  reflection  factor  for 
a  plane  wave  that  makes  the  name  angle  with  the  impedance  surface  of  the 
wedge  as  the  exit  angle  (p  . 
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For  negative  cp,  ty(-  ;  q>  +  rt)/t(cpQ)  =  ty(:fP!  -  n)/t(-Cf<o)  is  no  longer 
equivalent  to  the  reflection  factor,  but  is  similar  to  the  first  factor 
\|/(<i)  -  jt )  for  positive  cp.  At  the  shadow  boundary  of  the  incident  wave, 

<p  =  - (rt  -  cpQ),  cp  +  jr  =  cpo  and  ^ (cp  +  n )/i|/ (cpQ )  =  t(cpo)/\|/(cpQ)  becomes  equal 
to  unity.  The  first  factor 


i|/  ( -  I  cp  I  -  it ' 

V(fP  ) 


i(9’) 

-R-(cpoj  t{/(cpo) 


(8.11) 


where  R  (cp^)  is  the  exit  factor  for  the  second  surface  ('3  =  i3_),  then 
represents  the  reflection  from  the  second  surface. 

It  is  convenient  to  call  the  factor  the  shadow  boundary  factor. 
It  reduces  to  R(cpQ)  at  the  shadow  boundary  of  the  reflected  wave,  and  to 
unity  a.t  the  shadow  boundary  of  the  incident  wave , 
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IX.  ZEROES  AND  POLES  OF  THE  ^  FUNCTIONS  AND  THE  VARIATION  OF  THE  \|/ 

FACTORS  IN  THE  MALYUZHINETS  SOLUTION 

The  properties  of  the  wedge  surface  are  represented  by  the  ij functions. 
The  zeroes  and  poles  of  the  functions  follow  from  the  mathematical  deriva¬ 
tion.  They  can,  however,  be  derived  directly  on  the  basis  of  the  physics  of 
the  problem.  Since  the  Brewster  angle  adds  additively  to  the  argument  of  the 
functions,  we  can  assume  the  simplest  i-jossible  situation,  e.g.  matching 
for  normal  incidence  for  the  front  surface  impedance  (0+  =  0,  <Pq  =  <t>  -  jt  /  2  ) , 
and  pressure  release  behavior  for  the  rear  surface.  The  pressure  release 
part  then  is  very  large  and  independent  of  cp  and  <p  and  cancels  out  in  the 
ratio  '!'(Q')yV(cp  ) .  Thus  we  have 


b(q)  ^  y^-1-  - *>  _  (a) 

"Tr”  7/7v 


(9-D 


The  reflected  wave  in  the  solution  is  proportional  to  ty(jt  +  cp).  Since 
reflection  must  vanish  for  normal  incidence  when  <p  =  cp  =  <I>  -  jc/2,  we 
must  have 


^,/o(T  +  «)  1 


4‘/2  ««d  =  a-  -  n/2  ~ 


=  lU/J*  +  (*  "  «/2)]  *  0 


(9.2) 


Hence  the  zeroes  closest  to  a  =  0  cf  the  functions  ^(a)  arc  given  by 

+  ao  -  rt/2  +  21  (9.5) 

The  wedge  solution  contains  the  factor  cos  ( jfp  /2<J> ) ,  which  vanishes  for 
grazing  incidence  unless  ^  part  becomes  infinite.  For  a  rigid  surface, 

=  jt/2,  and  the  reflected  wave  never  vanishes,  not  even  at  grazing 
incidence  when  cp,  -■=  -  1  and  cp  =  <K  Time  (;(  H  1 )  must  become  infinite. 


kl 


Thus  we  have  with  G  =  jc/2 

v(a)-*  ^(a  +  ®  +  n/2)t(Ii(a  -  <f  -  jc/2)  » 


hence  at  cp  =  4> 


♦  («  +  « )— ^  ^(2®  +  ^y)  ^(jt/2)^,  00 


and  the  poles  must  be  given  by 


(9-M 

(9-5) 


+  a  =  2$  +  3it/2 


(9-6) 


The  ijr  functions  have  poles  and  zeroes  also  for  |a|  >  3«/2  +  20  ,  i.e. 
outside  the  basic  a  -  range.  But  these  poles  and  zeroes  are  of  interest 
only  for  the  product  and  series  developments  of  the  \|/  functions  (see,  for 
instance,  Eq.  A31). 

The  terms  ty  (cp  -  jc)  and  \|f  (cp  +  n )  of  the  stationary  phase  solution  are 
made  up  of  the  product  of  four  functions.  Each  term  is  zero,  when  one 


of  tuC-  arguments  of  the  iji  f UuC  tiuns  1-eduCes  to  ct  = 


.  r  ^  n 

■+■  14V  -r  — j 


9j  - 


xne  zeroes 


of  the  ^(ic  +  cp)  and  tU  -  cp)  terms  of  the  Malyuzhinets  solution  then  depend 

on  the  0  and  0_  of  the  imijedance  of  the  wedge  surfaces.  Let  us  first 
+ 

consider  the  term  \|/(n  +  cp)  which  in  the  space  0  <  cp  <  $>  is  similar  to  the 
reflection  factor  of  the  surface  z+.  The  argument  of  the  factors 


^(cn  +  it  +  <t  +  0+)  =  1)^  [(cp 


+  ~  -  O'  +  0+)  +  (24.  +  §)] 


(9-7) 


reduce  to  a  =  24>  +  ix/p,  ,  and  consequently  vanishes,  iff 


(9.8) 


i.e.  when  the  exit  angle  is  equal  to  one  of  the  two  Brewster  angles  0^  of  an 
impedance  covered  infinite  surface. 


Similarly,  the  factors 


'Mcp  +  n  «  J>  +  0_)  =  'll » [cp  +  +  3n/2  -  (?A‘  +  it/2)  +  G_] 


“  ^[9  +  |  -  30  +  (20  +  |)  +  OJ 


or  due  to  the  periodicity  of  t 


=  ^[cp  +  §  +  ®  +  6_  +(2$  +  |)] 


vanishes  iff 


9  =  -<!>-  3tc/2  +  G_ 


■which  is  not  in  physical  space,  or  if 


-  ir/2  +  0  <  _  0 


(9.9) 


(9.10) 


The  last  value  lies  in  physical  space  if  Re(0_)  =  it/2,  i.e.  if  the  z_  surface 
is  rigid  and  the  wedge  is  degenerated  to  a  semi-infinite  plane. 

The  first  term  in  the  Malyuzhinets  solution  is  proportional  to  \]c  (cp  -  n ) , 


The  zero  values  are  given  by 


11(1  -  it  a  <J>  +  0+)  -  1)1. [cp  -  it  -  <I>  +  (?A>  +  it/2)  -  it/2  +  GJ 


(9.11) 


Thus,  due  to  the  U<t>  periodicity  of  i|/  , 


cp  =  3n/2  +  4  +  0  -  Inf 


-  9  <  5<P  +  0  -  ~  <  <3. 

~  Q 


Thus  zeroes  in  the  first  factor  can  occur  only  if  <1  <  3it/U  +  0  /2.  It  is 
convenient  to  consider  the  space  qj  >  0  only  and  to  define  for  this  sj^ace 
the  factor  that  proportional  to  i|i(cp  1  it)  as  the  generalized  shadow 
boundary  (and  reflection)  factor,  and  the  factor  that  is  propiortional  to 
t(cp  -  k)  and  which  is  usually  much  smaller  than  the  scattering  factor. 
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The  term  ty(cp  -  it)  of  the  solution  then  is  practically  independent 
of  the  impedance  covers  of  the  wedge  as  long  as  cp  >  0.  We  have 

V(cp  -  tc)  =  ^(cp  -  it  +  $  +  0+)  -  it  +  <J>  -  0+)  ^(cp  -  it  -  d>  +  6J 

^(cp  ~  if  -  ®  -  0_)  ;  (9-13) 

all  the  factors  have  arguments  that  are  within  the  limits 

2$  +  Q+  -  it  «  0  +  0+  (9.H) 

and 

_  it  -  $  _  9_  &  -  2$  -  e_  (9.15) 

Their  product  therefore  is  practically  independent  of  0  ,  and  depends  on 
0_  only  as  cp  approaches  zero.  The  term  i|t(cp  -f  it)  is  taken  as 

ip(q>  +  it)  =  ^((p  +.  it  +  4>  +  ©  )  ^(cp  +  ir  +  4  -  0+) 

’  -  it  ~  0  +  0_ )  (cp  +  it  -  i>  -  e_)  (9.16) 

The  first  two  factors  depend  greatly  on  0+  because  their  zeroes  are  in  the 
range  of  positive  cp.  The  arguments  of  the  third  and  fourth  terms  then  are 
relatively  small  so  that  their  product  differs  only  little  from  unity.  Thus 
for  positive  cp,  ^(rp  -  u)  is  independent  of  and  depends  on  only  in  the 

neighborhood  of  cp  =  0.  In  contrast,  V(cp  H-  it)  depends  strongly  on  0+  and 
is  practically  independent  of  0  ,  For  negative  values  of  cp  =  -  |rp| 


\K(cp  +  «)  =  \J f(-  j  cp  I  h-  Jt)  =  4f(|fp  i  -  It) 

t(cp  -  n)  =  t(-  |  cp  {  -  n)  =  U/ ( j q> !  +  n)  (9-17) 


Changing  the  sign  of  cp  interchanges  the  two  terms  of  the  solution,  and 
the  second,  term  in  the  space,  cp  <  0  becomes  similar  to  the  first  term  in  the 
space  cp  >  0  except  that  0+  and  9_  are  interchanged.  But,  we  have,  however 
to  keep  in  mind  that  these  two  terms  are  not  independent  of  each  other  and 
that  strictly  speaking,  we  cannot  seperate  the  solution  into  reflection  and 
scattering  effects.  Figure  9.1  shows  the  two  ^  functions  as  a  function 
of  the  scattering  angle .  These  curves  have  been  computed  on  the  assumption 
that  the  two  surfaces  of  the  wedge  have  the  same  acoustic  impedance.  The 
functions  =  1)/ (cp  -  rt)  and  tyg  =  i)/(cp  +  it)  then  are  completely  symmetric 
with  respect  to  the  angle  cp  —  0,  i.e.  '('(-cp  -  ji)  =  vji(ip  +  ji),  and  'p(cp  -  jt)  = 
tc ( — cp  +  jt).  The  term  zero  ^or  Brewster  angles  of  exit  for 

cp  >  0  and  the  scatter  term  is  zero  for  the  Brewster  angles  of  exit 

with  respect  to  the  second  wedge  surface  cp  =  -  <J>  when  cp  <  0  . 


Figure  (9.2  )shows  the  i|/  functions  for  a  straight  edge  whose  faces 
have  different  acoustic  impedances.  In  Fig.  9.2a,  the  straight  edge  is 
matched  to  the  medium  for  normal  incidence.  The  shadow  boundary  factor 
ty(cp  -1  jt )  then  is  very  similar  to  that  for  a  straight  edge.  The  scattering 
factor  is  similar,  too,  on  the  side  of  the  incident  wave.  In  the  half  space 
opposite  to  that  of  the  incident  wave,  the  scattering  factor  is  modulated 
by  the  exit  factor  of  the  impedance  that  faces  the  shadow  space.  Fig.  (9.3) 
shows  the  same  as  some  of  the  ij/ (cp  T  jt )  and  \|/(cp  -  it)  functions  as  polar 
graphs . 
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X.  THE  EFFECT  OF  THE  IMPEDANCE  OF  THE  WEDGE  SURFACES 
ON  THE  DIFFRACTION  FIELD 
(a)  The  gackscattered  Field 

The  analysis  of  backscatter  is  simpler  than  that  of  diffraction  since 
we  restrict  ourselves  to  the  half  space  in  front  of  the  scatterer.  The 
dorninc nt  term,  also  in  the  backscattered  field,  is  the  second  term  of  the 
solution.  It  is  obtained  by  multiplying  the  angle  factor  Mp  with  the 
Malyuzhinetz  function  i|r(cp  +  it  )/i|/ (cpQ ) .  This  function  reduces  to  the  reflection 
factor  for  normal  incidence  and,  except  for  angles  near  cp  =  0,  approaches 
closely  the  exit  factor  of  the  reflecting  surface.  It  is  practically  in¬ 
dependent  of  the  material  on  the  other  side  of  the  straight  edge.  The 
backscattered  +  cpQ  )/iJ/ (fP0 )  field  is  thus  modulated  by  the  exit  factor 

of  the  wedge  material ;  it  vanishes  for  the  two  Brewster  angles.  But  in 
addition  to  this  field,  there  is  the  field  described  by  the  first  term  of 
the  solution  -  it)  [see  Eq.  (9-11)  )  and  Figs.  10.1  and  10.2).  In 

contrast  to  the  field  +  cP0)/'K('f)0),  this  field  does  not  depend  on  the 

impedance  of  the  illuminated  side  of  the  straight  edge,  and  depends  on  the 

impe  iance  of  the  face  in  the  shadow  side  only  for  small  angles  cp^.  We  would 

feel  tempted  to  interpret  the  Mp  term  as  the  field  generated  by  the  shadow 
boundary  of  the  reflected  wave,  and  the  M  term  as  that  due  to  the  shadow 
boundary  of  the  incident  wave.  However,  this  interpretation  of  the  Malyuzhinets 
solution  is  not  permissible. 

For  backscatter  in  the  direction  cp  --  0,  the  two  Malyuzhinets  components 
always  have  the  same  sign  and  add  up;  this  follows  by  considering  the  sign  of 
the  components  M^Kcp  -  «)  and  (T  +  n)  in  Figs.  10.1  and  10.2,  Thus  there 
is  always  a  maximum  of  backscatter  for  cpQ  -  0  like  in  the  case  of  the  ideally 
pressure  release  straight  edge  or  wedge  (see  Fig.  10,2).  The  width  of  this 
maximum  decreases  as  the  surface  impedances  increase  for  a  close  to  rigid 
surface;  all  that  is  left  is  a  needle  shaped  loop  and  for  a  rigid  straight 


edge,  the  width  of  this  needle  decreases  to  zero.  The  rigid  edge  does  not 
backscatter  for  cp  =  0.  The  term  (*P  T  it)  has  zeroes  at  the  two  Brewster 
angles  and  changes  its  sign  between  thin  lines.  The  term  M-hlf(cpo  -  jt)  is 
always  positive  from  cpQ  =  0  up  to  the  first  Brewster  angle.  Both  terms  then 
are  positive  and  the  field  is  represented  by  their  sum.  Beyond  the  first 
Brewster  angle  up  to  90° ,  the  second  term  becomes  negative  and  stays  negative 
until  cp  reaches  90°,  In  this  angular  range,  the  two  terms  substract.  At 
the  first  Brewster  angle,  the  resultant  field  is  finite.  But  above  the  first 
Brewster  angle,  the  two  terms  eventually  become  equal  in  magnitude  and  the 
resultant  field  is  zero.  A  computation  based  on  the  formulae  given  in 
Appendix  A  shows  that  this  zero  for  a  straight  edge  is  given  by  cos  0  +  cos  cpQ 
In  the  example  shown  in  Fig.  10.3,  the  zero  in  the  Brewster  angle  is  shifted 
by  about  27  degrees  towards  the  normal  of  the  surface.  At  and  above  the 
second  Brews  Lei-  angle,  both  terms  are  again  of  the  same  sign  and  their  cum 
decreases  drastically  towards  grazing  exit;  the  second  Brewster  angle  there¬ 
fore  does  not  seem  to  generate  another  zero  in  the  resultant  solution, 

(b )  The  Diffracted  Field 

In  considering  diffraction,  we  also  have  to  consider  the  solution  for 
negative  values  of  cp.  In  the  half  space  behind  the  diffracting  plane,  at 
the  shadow  boundary  of  the  incident  wave,  the  i)/(rp  4-  j()/\J/(cp  )  term  then 
reduces  to  unity  since  the  field  near  the  shadow  boundary  of  the  incident 
wave  is  not  influenced  by  the  material  of  the  reflector.  Away  from  the 
shadow  boundary  towards  the  positive  half  space,  this  field  shows  a  very 
slight  dependence  on  fl+.  In  contrast,  the  ^[^(cp  -  jt )/t|/ (rp^ ) ]  term  turns 
out  to  be  modulated  by  the  exit  factor  of  the  0_  surface  (see  Fig.  9*3-)  • 

Thus,  the  term  which  at  first  sight  we  would  have  expected  to  be  generated 
by  the  shadow  boundary  of  the  reflected  wave  is  modulated  not  with  the 


reflection  factor  of  the  reflecting  surface,  but  with  the  exit  factor  of 
the  surface  in  the  shadow  region.  We  are  compelled  to  assume  that  this 
term  is  generated  exclusively  by  the  surface  in  the  shadow  space. 

We  thus  arrive  at  the  conclusion  that  the  diffraction  field  of  a  rigid 
or  pressure  release  diffraetor  can  be  basically  attributed  to  the  shadow 
boundaries.  But  the  diffraction  field  of  an  impedance  covered  surface  is 
no  longer  generated  by  the  shadow  boundaries  alone,  but  is  strongly  modulated 
with  the  exit  factors  of  the  diffracting  surfaces.  It  exhibits  minima  in 
the  direction  of  the  Brewster  angles.  Thus,  the  orinetation  of  the  diffract¬ 
ing  surface  is  a  significant  factor  in  the  exact  theory  in  addition  to  the 
location  of  the  shadow  boundaries. 
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XI.  COMPARISON  OF  THE  MALYUZHINETZ  THEORY  WITH  THE 


KIRCHHOFF  AND  SOMMERFELD  THEORY 


In  the  calssical  Kirchhoff  theory  of  diffraction,  the  diffracting  screens 
are  black  screens.  The  screens  act  as  if  they  were  the  openings  of  infinitely 
long  ducts.  Any  radiation  that  hits  the  screen  disappears  in  it  without  the 
slightest  reflection.  Radiation  that  propagates  grazing  to  a  screen  is 
neither  attenuated  nor  diffracted  into  the  screen;  and  radiation  that  has 
been  diffracted  once  is  no  longer  affected  by  the  screen.  The  field  at  the 
aperture  is  assumed  to  be  that  of  the  incident  wave  and  diffracted  fields 
that  propagate  across  the  aperture  are  neglected  in  setting  up  the  solution, 
as  a  consequence  of  these  assumptions.  The  diffracted  field  spreads  mostly 
in  the  direction  of  the  incident  wave,  and  there  is  no  backscatter  in  this 
direction.  This  conclusion  follows  also  directly  from  the  classical 
Kirchhoff  integral.  This  integral  on  the  velocity  potential  U  contains  a 
source  and  a  dipole  term  as  follows : 


'  JhiiT  e'^  5u  1 


-  =  -  ~  (pcj  kU  cos  0  -  pcvn) 


1  e"^ 

—  (p  cos  9  -  pcv  )  ~i - 

pc  n  4«r 


(U.l) 


where 
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P  =  p  ZZ 


-  -dU 

n  bn 
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and  the  term  proportional  to  1/r  has  been  neglected  here.  For  a  plane  wave  such 
as  the  field  that  is  assumed  in  the  aperture  in  the  Kirchhoff  theory,  these 
two  terms  can  be  combined  to  a  term  of  cardioid  characteristics.  Because 


for  a  plane  wave  pcvn  =  p,  the  Kirchhoff  integrand  reduces  to 


!<9 


e-dkr 

-  p  (cos  0-1)  (11.3) 

Thus  the  radiation  is  zero  in  the  backward  direction  (0  --  0). 

The  diffraction  field  in  the  Kirchhoff  theory  can  be  attributed  to  the 
discontinuity  of  the  field  at  the  shadow  boundary.  This  shadow  boundary  acts 
like  a  semi-infinite  plate  that  vibrates  infinitely  slightly  below  its  coin¬ 
cidence  frequency.  The  regions  of  positive  and  negative  volume  flow 
counteract  each  other  except  for  a  region  of  about  half  a  wave  length  width 
along  the  aperture.  This  region  then  generates  the  so  called  edge  wave. 

Edge  radiation  and  radiation  from  the  shadow  boundary  are  thus  equinatural 
interpretations  of  the  diffraction  that  is  generated  by  a  Kirchhoff 
diffractor.  it  turns  out  that  for  a  straight  edge  the  Kirchhoff  diffraction 
is  a  maximum  in  the  direction  of  propagation  of  the  incident  wave,  and 
exactly  zero  opposite.  It  is  crudely  inversely  proportional  to  the  angular 
distance  from  the  shadow  boundary  and  hence  decreases  monotonically . 

The  Sommer f eld  edge-wedge  theory  applies  for  rigid  or  pressure  release 
straight  edges  and  wedges.  The  first  term  of  the  Sommerfeld  solution 


u(r,<$) 


1 


cos 


(11.M 


is  independent  on  whether  the  diffracting  surface  is  rigid  or  is  a  pressure 
release  surface.  It  satisfies  the  wave  equation  and  satisfies  exactly  the 
same  conditions  at  the  screens  as  the  Kirchhoff  solution.  In  the  Sommerfeld 
language,  the  screens  represent  the  branch  surfaces  (the  two-dimensional 
branch  cuts)  that  connect  the  physical  sj>aee  with  the  Riemann  spaces  (the 
ducts  that  represent  the  Kirchhoff  screens).  It  satisfies  the  continuity 
conditions  at  the  shadow  boundary  of  the  incident  wave.  But  the  Sommerfeld 
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term  does  represent  the  exact  field  in  the  apertures  as  it  would  correspond 
to  Kirehhoff  screens  if  they  could  be  realized.  Because  of  this  field, 
backscatter  is  quite  considerable. 

The  corresponding  Malyuzhinets  solution  is  obtained  by  assuming  the 
impedance  matched  for  the  angle  of  incidence  cpQ  =  «/ 2: 

cos  ©  -  cos  0 

£  _  Q  _  O _ 

COS  ©  +  cos  0 


(.11.5) 


3C 

Hence  0  =  +  cp  =  —  ;  the  negative  sign  is  unacceptable  because  the  impedance 

o  ^ 

must  have  a  positive  real  part.  This  angle  cpQ  'has  been  introduced  in 
Section  5  as  the  acoustic  Brewster  angle. 

Figure  11.1  and  11.2  show  a  comparison  for  normal  incidence  of  the  Kirehhoff 
solution  with  the  first  term  of  the  Sommerfeld  solution  and  with  the 
Malyuzhinets  solution  for  impedance  match.  The  Kirehhoff  solution  shows 

no  backscatter,  but  otherwise  agrees  better  with  the  exact  solution  than 
the  first  term  of  the  Sommerfeld  solution.  The  zero  of  the  Malyuzhinets 

solution  at  rp  =  n,  i.e.  at  the  impedance  matched  surface  is  obvious.  .Since 
we  have  used  a  stationary  phase  evaluation,  the  distance  from  the  edge  of 
the  wedge  is  assumed  to  be  very  large,  as  a  consequence,  all  the  field 
energy  is  impedance  matched  at  and  near  the  impedance  matched  surfaces.  For 
finite  wedges  and  edges,  the  zero  will  usually  be  replaced  by  a  minimum.  It 
is  very  likely  that  the  Kirehhoff  solution  then  will  lead  to  still  a  much  better 
agreement  with  the  exact  solution  that  it  does  for  the  infinitely  extended  wedge. 

Figure  11.3  shows  the  Wo  terms  of  the  Malyuzhinets  solution  plotted 
separately  and  their  sum,  for  impedance  match  at  normal  incidence.  Tlie 
first  term  1  is  negative  and  reduces  the  resultant  amplitude  to  zero  for 
rays  that  prorogate  grazing  to  the  impedance  matched  surface. 
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In  the  Kirchhoff  theory,  diffraction  is  exclusively  a  phenomenon  of 
the  shadow  boundaries  or  waves.  We  therefore  obtain  the  Kirchhoff  solution 
by  simply  adding  the  field  generated  by  the  incident  wave  and.  that  generated 
by  the  reflected  wave,  assuming  that  the  reflected  amplitude  is  R(<p  )  times 
that  of  the  incident  wave.  This  solution  then  approaches  closely  the  exact 
Malyuzhinets  solution  on  the  whole  angular  range  except  for  two  regions : 

(l)  the  Malyuzhinets  solution  vanishes  along  the  impedance  covered  surface 
because  we  have  derived  it  by  a  stationary  phase  method,  assuming  that  the 
field  point  is  very  far  away  from  the  diffracting  edge.  In  contrast,  the 
Kirchhoff  surfaces  do  not  affect  grazing  rays .  The  Kirchhoff  and  the  exact 
solution  therefore  must  differ  for  cp  =  +  $;  and  (2)  the  Kirchhoff  terms  always 
lead  to  zero  baekscatter.  The  two  solutions  with,  therefore,  also  differ 

for  q>  -  -r  r.i  . 

O 

The  Malyuzhinets  solution  shov/s  that  the  exit  factor  should  be  used 
rather  than  the  reflexion  factor.  For  the  exit  angle  equal  to  the  angle 
of  incidence,  the  two  factors  arc  the  same.  But  for  other  angles,  the 
results  differ  slightly.  The  exit  factor  usually  leads  to  a  better 
approximation  of  the  position  of  the  minima  of  the  diffracted  field;  but 
in  general,  the  deviation:;  are  not  great  because,  also  in  the  view  of  the 
exact  theories,  diffraction  is  basically  a  phenomenon  of  the  shadow 
boundaries.  It  is  the  angle  of  incidence  and  that  of  reflection  that 
determines  the  bulk  of  the  diffraction  phenomenon.  Figures  11  Ji ,  11  .<)  and 
11  .6  show  a  comparison  of  the  exact  solution  with  the  Kirchhoff  and  Gomrnur- 
feld  solution,  the  second  terms  weighted  both  with  R(<j>  )  and  ]< (cp).  Agree¬ 
ment  between  the  exact  solution  and  the  Kirchhoff  result  is  excellent  for 
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the  regions  pointed  out  above  i,e.  grazing  exit,  backscatter  in  direction  of 
incidence  and  reflected  waves. 

Tinally,  we  are  very  interested  in  wedges  where  one  surface  is  rigid, 
the  other  pressure  release.  The  ^  function  for  this  case  reduces  to 


♦  (c)  =  +  «  +  f)+4(a  ♦  ♦-§)  =  ♦!(§)  eos(a  +  ♦)  &  =  +|(f)  «s(g  +  J) 


(11.6) 


where  we  have  assumed  the  +  $  surface  rigid,  the  -  $  surface  pressure  release. 
The  functions  (a  +  )  have  been  dropped  since  they  are  constant  and  occur 

V 

in  similar  form  in  ty(rp  ).  The  results  are  shown  in  Figure  11. 7.  If  the 
wave  impinges  on  the  pressure  release  surface,  and  if  the  second  surface 
is  rigid,  the  solution  is  similar  except  that  the  term  n/k  is  replaced  by 
-  jt/'i .  The  result  of  the  computation  is  also  plotted  in  Figure  11. 7. 

For  normal  incidence,  the  curves  are  practically  independent  of  the 
impedance  of  the  material  in  the  shadow  space.  It  follows  from  the  derivaliion 
in  section  (  IX ),  that  the  solution  near  the  shadow  boundary  is  entirely 
independent  of  the  surface  properties,  and  that  it  does  slightly  dej>end  on 
the  impedance  of  the  surface  in  the  shadow  space  as  the  field  point  moves 
towards  the  diffracting  plane.  As  a  consequence  of  this  dependency,  the 
zero  in  the  diffracted  field  which  for  Cj;  =  jc/2  occurs  at  exactly  zero  degrees 
for  the  rigid-rigid  straight  edge  is  shifted  to  about  35°  for  the  rigid- 
pressure  release  straight  edge.  Figure  11.8a  shows  the  exact  solution  for 
the  rigid-rigid  straight  edge.  Fig.  11.8b  that  for  the  pressure  release- 
release  straight  edge.  Figure  11. y  shows  for  comparison  the  Kirehhoff 
approximation,  which  as  would  be  expected,  leads  to  very  good  agreement  in 
most  of  the  angular  range. 


Figure  11.10  shows  the  exact  solution  for  a  9°°  pressure  release- 
rigid  rigid,  and  a  rigid-pressure  release  wedge.  For  comparison,  the 
diffraction  at  similar  rigid-rigid  and  pressure  release  pressure- 
release  wedges  are  also  entered.  Again,  the  Kirchhoff  theory  would  lead 
to  a  very  good  approximation  of  the  exact  curves.  For  instance,  if  the 
incident  wave  impinges  on  the  pressure  release  surface  at  cp  =  120°,  its 
shadow  boundary  occurs  at  cp  =  -  6o°.  The  shadow  boundary  of  the  reflected 
wave  which  now  is  reflected  in  antiphase,  occurs  at  cp  =  -  30°,  The  results 
for  an  incident  angle  of  135°  are  not  included.  This  position  lies  along 
the  face  of  the  wedge.  The  solution  for  this  case  does  not  exist.  As 
stated  previously,  the  stationary  phase  method  assumes  that  the  field  point 
is  very  far  from  the  diffracting  edge,  and  since  the  angle  of  incidence 
lies  along  the  surface  of  the  wedge,  all  the  incident  energy  is  absorbent 
for  these  rays  propagating  at  grazing. 
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XII.  THE  PHYSICAL  INTERPRETATION  OF  THE  HALYUZHINETS  SOLUTION 

The  exact  solution  for  a  straight  edge  with  non  rigid  surfaces  is  based 
on  the  solution  for  the  straight  edge  with  pressure  release  surfaces.  There 
is  a  good  reason  for  using  the  pressure  release  solution  as  a  starting  point 
and  to  modify  this  solution  for  impedance  covered  surfaces.  A  wave  pro¬ 
pagating  in  the  direction  of  the  plane  of  the  edge  towards  the  edge  of  a 
pressure  release  straight  edge  is  backscattered,  i.e.  the  diffracted  field 
is  not  zero  for  cpQ  =  0.  It  is  physically  obvious  that  an  impedance  covered 
straight  edge  will  also  backscatter,  whereas  a  rigid  straight  edge  will 
have  no  effect  at  all  on  a  wave  that  propagates  parallel  to  its  plane, 

being  zero.  The  wave  passes  left  and  right  along  it  without  generating 
a  diffraction  field.  Thus  a  solution  based  on  the  field  scattered  by  a 
rigid  straight  edge  would  break  down  for  cp  =  u.  The  so3.ution  is  finite  for 
cp  =  +  $  for  the  pressure  release  wedge,  and  for  all  wedges  whose  faces  have 
finite  impedances  regardless  of  how  small  their  impedance  may  be.  The 
vanishing  of  the  basic  solution  for  cp  =  +  <v  is  therefore  of  no  consequence. 

The  solution  for  the  pressure  release  wedge  then  is  given  by  the  sum 
of  the  two  Sommer f eld  terms,  or  by  the  sum  of  the  two  Malyuzhinets  terms. 

The  Sornmcrfeld  and  the  Malyuzhinets  terms  represent  different  partial 
fraction  representation  of  the  exact  solution  for  the  pressure  release 
straight  edge.  The  Sornmcrfeld  solution  is  equivalent  to  superposing  the 
fields  generated  by  a  source  and  by  an  image  source  located  behind  the  com¬ 
plementary  straight  edge. 

However,  the  image  concept  breaks  down  when  the  surfaces  are  impedance 

2  3 

covered.  An  infinite  number  of  image  sources  of  strength  R,  R  ,  R  ,  etc. 
would  then  be  needed  where  R  is  the  reflection  factor  of  each  of  the  two 


surfaces,  and  the  Somrnerfelu  representation  becomes  useless. 
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The  Malyuzhinets  term  M^1!'  (o  +  jt)/^(cpo)  depends  almost  only  on  Q+; 
and  it  depends  only  to  a  small  extent  on  the  impedance  of  the  surface  on 
the  shadow  side  of  the  straight  edge*  in  the  angular  range  near  zero.  It 
reduces  to  the  amplitude  of  the  reflected  wave  at  the  shadow  boundary  of 
this  wave,  and  to  unity  at  the  shadow  boundary  of  the  incident  wave.  For 
other  angles  on  the  positive  side  of  the  straight  edge,  this  field  is 
proportional  to  the  exit  factor  of  the  reflecting  (upper)  surface.  The 
field  described  by  this  term  is  not  very  different  from  the  Kirehhoff 
fields  generated  by  the  shadow  boundaries  of  the  reflected  and  the  incident 
wave.  But  as  we  penetrate  into  the  shadow  region,  the  contribution  of  the 
shadow  boundary  of  the  reflected  wave  practically  disappears. 

In  the  region  around  cp  =  0,  the  Malyuzhinets  term  ^(cp  +  « (<P0 )  describes 
a  field  that  is  stronger  than  the  Kirehhoff  field.  Because  of  the  impedance 
cover,  the  surfaces  of  the  straight  edge  move  in  and  out  like  the  surfaces 
of  a  transducer;  they  are  forced  to  vibrate  by  the  incident  wave.  If 
there  were  no  edge  region,  the  impedance  covered  surface  would  reernit  the 
incident  wave  in  form  of  the  reflected  wave  and  would  simply  act  like  a 
reflector.  However,  the  wave  that  passes  the  edge  rexiresents  a  reservoir 
of  energy.  Some  of  the  energy  is  also  sucked  into  the  edge  region,  a  wave 
is  always  bent  into  an  iinxieuance  surface,  and  remitted  as  an  edge  wave .  The 


^(rp  -  it)  =  ty  (cp  +  y+)irr(q>  -  -  2k  -l-  0_  )^(cp  -  2k  +  0  ) 

=  t  (p  -  2it  -1-  0_)i|/  (cp  -2 n  -  0  )  for  cp  >  0, 


and 


because  csJ  1  if  fa  I  <  n 
jc  •  — 


t|/(cp  -  it)  =  ty(-  |  cp  |  -  it)  "=  t  (  f  cp  |  +  2n  +  )t  (  j  cp  |  -t-  2 it  -  0  )  for  cp  <  0 
also 

^  (cp  +  it)  =  i|c  (fp  -(■  2  it  +  0+)i|/jt(rp  +  2  n  -  0+)+K(q»  d  0_ )  ^(qi  -i-  0_) 

"=  i}1  (fp  +  2it  +  0  )>p  (cp  +  2it  -  0  ) 

It  ’«■  3t  T" 
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first  term  of  the  Malyuzhinets  solution  then  seems  dominantly  to  represent 
two  fields.  The  function  \|r(cp  -  n)  depends  almost  exclusively  on  6  .  In 
the  front  space,  i|r(cp  -  jt)  depends  very  weekly  on  and  is  practically 
independent  of  0  .  This  field  can  therefore  be  interpreted  as  a  further 

t" 

contribution  of  the  shadow  boundary  of  the  incident  wave  to  the  diffraction 
field.  In  the  space  cp  <  0,  this  term  becomes  proportional  to  the  exit  factor 
of  the  surface  0_.  Hoe  term  H.  1  (cp  -  :t)/fy(cpo)  can  then  be  interpreted  as 
the  field  that  is  reemittcd  by  the  vibrating  edge  region.  Figure  (11. 5) 
shows  a  schematic  representation  of  the  nature  of  the  two  terms  of  the 
solutions  on  the  two  sides  of  the  diffracting  semi-infinite  plane. 

The  exact  solution  proves  that  the  Kirchhoff  concept  of  the  diffraction 
phenomenon  being  solely  generated  by  the  shadow  boundaries  is  only  a  first 
approximation  for  impedance  covered  surfaces.  According  to  the  exact  theory, 
some  of  the  incident  energy  is  scattered  at  the  edge  because  of  the  dis¬ 
continuity  of  the  acoustic  properties  of  the  medium  due  to  the  impedance 
cover.  Because  of  its  finite  impedance,  the  edge  region  is  compressible 
and  scatters  sound  like  a  chain  of  air  bubbles,  located  at  the  edge, 
scatter  sound  in  water.  This  scattering  effect  seems  to  be  greater  the 
smaller  the  surface  impedance .  A  pc  matched  edge  thus  scatters  almost  as 
much  sound  as  a  pressure  release  wedge;  for  instance,  when  the  wave  pro¬ 
pagates  into  the  edge  parallel  to  an  impedance  covered  plane.  In 
addition,  the  edge  region  and  the  whole  surface  of  the  straight  edge 
reemits  the  incident  energy  that  is  not  scattered  according  to  the  exit 
law.  The  amplitude  reemitted  ray  is  proportional  to  the  exit  factor  of 
the  surface.  The  diffracted  Brewster  angler,  of  the  impedance  covered 
surfaces,  regardless  of  the  position  of  the  shadow  boundaries.  If  the 
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Kirchhoff  concept  would,  apply,  then  the  diffraction  field  would  be  pro¬ 
portional  to  the  reflection  factor  for  the  angle  of  incidence  of  the 
radiation  and  would  not  show  any  Brewster  angle  dependence. 

The  curves  presented  in  this  report  show  the  variations  of  back- 
scatter  and  diffraction  in  those  instances  that  might  be  of  practical 
interest.  The  various  approximations  are  of  considerable  interest  since 
they  help  to  understand  the  physics  of  diffraction.  But  on  the  basis  of 
the  cosine  representation  of  the  Malyuzhinets  functions,  practically  exact 
computations  can  be  performed  with  the  same  case  as  Sommerfeld  or  Kirchhoff 
approximations.  There  is  thus  no  need  to  investigate  wedge  or  edge 
diffraction  with  the  aid  of  approximate  theories. 
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APPENDIX  A 

Representation  of  the  i|/$  Functions 

The  following  representations  of  the  \)^(a)  functions  have  been  derived 
in  the  above  Section  V. 


The  form  A1  is  derived  in  very  much  the  same  manner  as  Eq.(5.l8)  except  that 
the  w  integration  is  performed  first  in  Equation  (A2).  By  deforming  the  path 
of  integration  from  infinity  above  the  real  axis  to  zero  and  back  to  infinity 
below  the  real  axis,  and  determining  the  residues,  the  following  forms  aye 

obtained  for  the  rational  ratio  h<b/?.  =  n/m  and  a(k,f)  =  rr  (— - — — — ) 

we  have  for  odd  and  even  n,  respectively: 


/  v\  _  n  77  1  cos  1/2  afk,,1?)  N 

^in/f4m  cos  1/2  [a/n  +  a(k)j)]'/ 


(-1) 


jin/tiri 


\n/^(a) 


k-1  i- 1 


n  l\  7  xa(k,iH  a/n 

"  “»L<- '  y7 

a(k,i) 


u  ct-g  u  du 
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Examples : 


V(a)  =  cos  i ;  V(a)  ~  exp[^t/'u  ~issu?“ du]  •  (A6) 

O 

^  (a)  =  COS  1/6  (g  -  rt)  COS  1/6  (a  4-  n)  .  (A7) 

cos2(jt/6)  cos  1/6  a 


^(a) 


=  exp 


-  a 

if 


-it  sin  v  +  cos  £■ 


Sin  2 


-  2v 


COS  V 


dv 


•# 

m 


(AS) 


For  large  values  of  Im(z): 


HZ 


exp 


2  it 


/ 


in  cosh  ( jts /U<v ) 


cosh  s 


d.3  ] 


(A9) 


and 


i teL  ->  i 

*<V 


(A10) 


The  following  relations  can  be  derived  for  the  Malyuzhinetz  functions  either 
from  their  integral  representation  or  residue  evaluation  of  the  contour 
integral  or  from  the  basic  properties  of  these  functions  (moromorphic  in 
strip,  poles  and  zeroes  known). 


where  z 


6o 


/  jt  \  f  it \  1 2  / Tt \  JtQ> 

V°  +  Z]  f$(a  ‘  2}  =  %  (2}  C°S  E* 

(All) 

va)  va  +  rt)  =  4(f) cos  & (a  +  f) 

(A12) 

Va  +  ^  Va  -  =  V2(a) 

(A13 ) 

+4(a  +  2®)  ^(a)  =  ^(S)2  ^a,/2^a  +  ^ 

* 

(Al4 ) 

Va  +  2$)  ]_  * 

2»J  =  «*  2  (“  +  2>  > 

(A15) 

V  (a  +  4®)  l  ,  t 

'  »t(a)~  »  =ot  f  (a  +  5  +  2»)  , 

(A16) 

V“'  j-j-r  =  tan  |  (a  +  |  -  2*) 

^(a  -  U<b)  2  2 

(A17) 

**(“  +  f )  Va  “  I?  =  4(f}  cos  §  ’ 

(A18) 

Va  +  *)  Va)  =  4(l}  cos  (a  + 

(A19) 

♦*(«)  va  - ,r)  =  4(f}  cos  A (a  -  f } 

(A20) 

1f9[z  +  (2®  +  —)]  -  +  sin  -Y-ro— ^  cocec  ^  (2®  -  g  +  ®) 

(A21 ) 

ap  =  ±  t24,  H  f1)  >  zo  =  i  (2a(  H'  f ] 

(A22 ) 

the  values  of  z  for  the  poles  of  ^(z) 

and  z  those  for 
0 

the 

zeroes . 
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The  residues  r,  at  the  poles  a,  are  easily  determined  by  means  of  the 

+  "T 

relation 

Vz  t  (2i>  +  T)]  =  i  sin  z)~ cosec  ff  V2$  -  §  ±  *)  >  <A23) 

which  puts  the  poles  into  a  trigonometric  factor, 

Tito  more  forms  for  the  complex  t(a)  function  are  of  interest.  They  are 
derived  with  the  aid  of  the  relations  that  apply  for  the  ^(a)  functions. 


t(a)  =  +  i  +  e+)  ^(a  +  ®  -  0+)  ^(a  -  +  e_)  ^(a  -  4>  -  6_) 

=  Cv#(|)]  [cos  5^  (a  +  *  -  ■fl+)  cos  (a  -  S>  +  *_)] 

(cc  +  -  —  jt  +  'fl  ")  * .  (a  -<!>  +  —  n  -  '5  ) 

x  - s - n  ou - 1 - -1  a ik) 

♦t(a  *  *  -  j  <  -  <t)  *4(a  -  is  +  i  ji  +  O 
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LL 


Hence 


l|/(cp  -  «)  =  COS  ~  (cp  -  It) 


i|f(cp  +  3t)  =  cos  ^  (cp  +  n) 


and  for  <I>  =  n: 


t(cp  -  jt)  =  const,  sin  cp/2 
\'c(q>  4-  jt)  =  const,  sin  cp/2 


If 


V  -  fl-  *  t + 


j/o) 

V(cpQ) 


=  1  regardless  of  a  and  cp 


If  one  surface  is  pressure  release 


6  =  f  f  >  td,(Q  + 


and  the  two  ^  factors  that  contain  0  drop  out  of  the  solution 


jK_cp+jO 

'KO 


If  a  surface  is  rigid,  0  -  jr/2  ,  then 


2  ,n  > 


V« +  * +  f)Va  +  *  -  !>  =  K  (t}  cos  ff 


(A2  6) 


(A27) 


o.nd  the  two  terms  in  the  solution  t(cp  +  jr)/'KcP0)  can  he  replaced  by  the  two 


terms  cos 


n  respectively. 
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1 


If  the  front  surface  is  impedance  matched  (z  =  pc),  then  0  =  0,  and: 


't'(a)  =  +  $>1 Iq(cl  -  4>)^(a  -  &  +  0_Ho(a  -  <2>  -  0_) 

=  ^(®)+4/2(a)'|fC»(a  "  4>  +  G-)^(a  -  C’  ~  O 


/•♦■pfc,  -.iM.fl  . 
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APPENDIX  B 

Comparison  of  the  Malyuzhinets  and  the  Gommerfeld  Solution 

The  Sommerfeld  solution  applies  only  to  the  rigid  and  the  pressure 
release  wedges.  It  consists  of  two  terms,  one  of  which  represents  the  field 
of  the  incident  wave,  the  other  the  field  generated  by  an  image  source  (hidden 
in  a  Riemann  space).  It  is  of  interest  to  force  the  Malyuzhinets  solution 
into  a  similar  form,  to  investigate  whether  in  the  general  case  a  simple  analysi 
is  possible  in  terms  of  fields  generated  by  a  real  source  and  an  image  source. 
The  Malyuzhinets  solution  is  given  by 


For  the  straight  edge,  v  =  2.  The  two  angular  functions  in  (BE)  can  he  written 
as  follows ; 


Let  us  write  their  sum  with  a  common  denominator. 
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For  the  minus  sign  we  have  the  numerator 


94-9  cp  -  9 

it  ,  o  it  ,  o 

cos  —  +  cos  -  -  cos  —  +  cos  - - - - 

V  V  V  V 


For  the  plus  sign  we  get 


9  +  9  9-9 

o  .  o 

COS  -  4  COS  - 


9 

o  ?  O 

2  cos  —  cos  — 
V  v 


(B5) 


COS  —  +  COS  - 

v  v 


9  +  9  9  -  9  , 

0  ,  it  o 

4-  cos  —  -  cos  - - 

V  V 


<J> 


=  2(cos  ~  -  sin  —  sin  2)  (b6) 


Thus  we  have 


£!  '  S2  = 


9 

O9o 
2  cos  —  cos  — 
V  V 


9-9  0  +  9 

/  ^  Ow  It  J  o 

(cos  —  -  cos - '(cos  —  4-  COS - 

V  V  V  V 


(B7) 


and 


S1  +  S2  = 


9 

„  /  11  ,  o  .  9  x 

2  (cos  —  -  sm  —  sin  —  J 
v  v  v  ' 


9-9  9  +  cp 

/it  Ox  /  rc  .  o 

(cos  —  -  cos  - ;(cos - 1-  cos - 

V  V  V  V 


(BO) 


The  Malyuzhinets  solution  contains  the  terms  —  and  —  where 


M 


1 


9  -  it  .  cp  4-  m  ^o 

sin - sin  —  sm  ' - -  sm  — 

v  v  v  v 


9  -  «  -  9, 


2  sin 


2V 


cot 


9, 


9  -  IC  4  9 

_ _ _C 

2v 


9  4-  It  -  9 


9  +  it  +  9 


2  sin 


cos 


o 


.  /Cp  +  TC  \  .  /Ox  r  CD  —  IT 

sin  - - )  -  sin  ( — )  -  sin  - 

v  '  v  v 


9.. 


im 


°-l 


9  -  it  -  9  v  +  rr  -  o  9  -  it  4  q>  c-  4-  n  4-  cp 

Y  .  0  o 

sm  - - - -  2  cos  - r — ■ — -  cos  - 


2  V 

(B9) 

(BIO) 


2  sin 


2v 


2v 


2v 


2V 


mtmsrnimumimumi  nmm/muma 


•  -.m 
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[cos 


[sin  -  "  ' -  sin  — ] 
_  v  v 


9 


CD 

r  9  ~  It  .  0"i 

Lsm - sxn  — J 


9  -  cp  cp  +  cp 

[cos  —  -  cos - — ]  [cos  —  +  cos - -■] 

V  V  V  V 


r  .  cp  +  it  ,  cp  _  it  . 

[sxn  __  _  sin  ___  ] 


„  .  it  cp 

2  sin  —  cos  — 

V  V 


9  +  cp 


cp  -  cp 

It  O-I  r  It  '  '  O I 

V  "  cos  T“J  [cos  V  H  cos  - - ] 


Similarly  we  have 


(cos  - —  cos 

V  V 


cp  -  cp  cp  +  cp 

Ow  It  ,  Ov 

-)(cos  -  +  cos  — - - ) 


(Bill 


JL  +  = 

Cl  6..  r?  -  it 


4- 


M  M  sin 


9_ 


-  sin  9  ,9+n  .  o 

v  o  sin  -  -  sm  — 


9 


[sin 


9 


•  w  i  I*  t  .  w  | 

sm  — J  4  Lsm  — - •  -  sm  —  J 

v  v  vJ 


9  -  9  „  94-9 

/It  0\/tt  o  \ 

(cos  —  -  cos - Jlcos  —  +  cos - ) 

'  -\J  -W  1 !  .1  ' 


9 

„  r  •  o  .  .  9  it  -i 

2  L-  sin - H  sm  —  cos 

v  v  vJ 


9-9  cp  -  9 

- -)(cOS  -  +  COS - ) 

V  V  V  / 


/  A 

(col;  -  -  cot 


(B12) 


We  can  er.press  the  Malyuzhinets  terms  in  tonns  of  the  Sommerfeld  terms  as 
follows : 


1_  _1_ 

f/’M 


„  .  Tt  CP 

2  sm  —  cos  — 
v  v 


(cos - COG 

V 


cp  -■  9 


)  (cog  ■—  +  COG 
V 


9  +  9, 


it  cp 

sm  —  cog  — 
v  v 


...  -t  VJ 

2  COG  —  cos  - 

V  V 


V 


:in  I 

~  W  ~  S2^ 


COG 


O 

V 


( 1.(13) 
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°M 


of  A  o  ^  .  <P  IT  o 

2  I  -  sin  —  +  sm  —  cos  — j 

_ _ _ v _  v  v _ 

cp  -  ffl  Cp  +  Cp 

,  jr  O  \,  «  ,  x  o, 

(cos  -  -  cos  — )  ( cos  -  +  cos  ^ ) 


/  •  O  ,  ,  Cp  Jtv 

(-  sm  —  +  sin  —  cos  — ) 

- - y - y — y-  (Sl  +  s  ) 

cp  _  v  1  2 

/  «  .  o  ,  cp  . 

(cos - sin  —  sin  —  ) 

v  v  v  v 


(EliO 


The  Malyuchinets  solution  can  now  be  written  as  follows 


i|;(cp  -  n)  +  'p(cp  +  it)  =  A 


(B15) 


ip(cp  <-  jt)  -  Ip(cp  +  «)  =  B 


(Bl6) 


■’K?  -  ip  _  »(<P  +  tc)i  _  A.  ,_1 _ L  B  ,1_  1_  v 

°M  2  °K  2  °M 


(El  7 ) 


C  9  (-  sin  —  +  Sin  ^  cos  £) 

u(r-P)  =  h^n\k  sin  v  (si  -  s2}  +  E  cos  ~ - - - —  (si H- s 

^  (cos - sin  —  sin  — ) 

v  v  v  (bi8) 


f  /  f  O  ■ .  Cp  JT  \ 

n  (  9  (-  sin - y  sm  -  cos  — ) 

—T—r <C(A  sin  —  +  B  sos  -2. - 1 - —S - 1-  ]  s 

v  v  *  .  .  <0  1 
^  cos  —  -  sin  —  sm  — 


+  t-A  sin  — t  — 

v 


cp  Cp  m 

t  0\t  .  O  .  Cp  Us 

v  ‘ '  v  v  v 


JT  .  O 

cos - sin  — 

v  v 


(B19) 


For  the  pressure  release  wed ge,  A/ip^  )  -  1,  B/ip(cl;o)  =  0,  and  the  solution 
becomes  identical  with  that  given  by  Sommer fcla.  For  the  straight  edge, 


this  Kq.  I9  reduces  to 


u(r/p)  = 


(  cos(cp  /2)  cos(9r/2)  *\ 

"  3  '.kM?  ^ 4  3  TEifm  1  y 
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it)  +  \|r(tp  +  Jt)  +  [t(9  -  Jt)-T|f(q>  +  It)] 


C°S  q)o/2\ 
“sin  q>j2J 


[_!j/(cp  -  it)  -  ijr(cp  +  tc)  +  [\|;(cp  -  jt) 


cos  cp  /2 

+  -STW21 


_ 1_ _ 

cp  +  cp 

o 

C°3  2 


(B20) 

The  first  factor 


(B21) 


has  a  pole  at  cp  =  -jt  +  cp^;  this  pole  occurs  at  the  shadow  boundary  of  the 
transmitted  wave.  The  second  factor 


has  a  pole  at  cp  =  jf  -  cpQ  at  the  shadow  boundary  of  the  reflected  wave. 

For  the  shadow  boundary  of  the  incident  wave,  cp  =  -  n  +  cpQ.  The  first 
bracket  with  the  factor  in  front  reduces  to  the  Sommerfeld  first  term:  the 
second  bracket  with  the  factor  in  front  reduces  to  the  shadow  boundary  of  the 
reflected  wave  to  the  classical  reflection  factor  (l/2  t(2jc  -  cp^V'HT^) 
times  the  second  Sommerfeld  term.  But  the  term  differs  from  the  Sommer le Id 
terms  in  the  range  between  the  two  shadow  boundaries.  This  all  means  that 
the  impedance  surfaces  generate  some  interaction  between  the*  two  fields  and 
that  the  exact  solution  therefore  cannot  be  obtained  by  the  method  cf  simple 
image  sources.  The  results  show  that  the  MalyuzhinoLs  form  of  the  solution 
is  considerably  simpler.  Forcing  it  into  the  Sommerfeld  form  helps  in  no  way. 
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Figure  11.  Sb  The  Sommerfeld  solution  for  diffraction  from  a  pressure  release  e 

identical  with  the  field  when  both  surfaces  are  pressure  release. 
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The  Kirchhoff  solution  for  diffraction  from  a  pressure  release  edge.  The 
rear  surface  does  not  contribute. 


Diffraction  from  a  90  wedge  for  d>  =  90  .  The  wedge  surfaces 
satisfy  the  various  combinations  of  rigid  and  pressure  release 
boundary  conditions  indicated  above. 
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G.D.  MALYUZHINETS 

Moscow,  Akusticheskiy  Zhurhal,  No  4,  1969.-  pp  629-630 

Excerpts:  Prof  Georgiy  Danilovich  Ma.lyuzhinetz,  Doctor  of 
Phys icomathematical  Sciences,  died  suddenly  on  14  August  1969 
at  the  age  of  59.  Malyuzhinets  was  a  prominent  Soviet  specialist 
m  acoustics  and  theory  of  wave  diffraction,  member  of  the  Communist 
Party  of  the  Soviet  Union,  and  head  of  the  laboratory  of  the  Acoustics 
Institute  of  the  USSR  Academy  of  Sciences.  .  .  . 

Malyuzhinets  was  very  acti've  in  dir  acting,  graduate  students  and 
in  giving  lectures  both  in  schools  (Moscow  State  University,  Moscow 
Physicotechrical  Institute  and  in  research  institutes.  His  comprehen¬ 
sive  approach  to  scientific,  problems  plus  his  great  intuition  in  physics 
attracted  the* -youth  to  him.  Many  young  scientists  did  their  disserta¬ 
tions  under  his  direction  or  with  his  constant  advice. 

Malyuzhinets  attached  exceptional  importance  to  the  training  of 
specialists  in  diffraction  theory.  lie  prepared  and  gave  an  original 
course  on  diffraction  theory  for  students  in  the  Moscow  Physico- 
tecbnioal  Institute.  He  strove  to  set  up  diffraction  theory  as  a 
speciality  in  the  Mechanics  and  Mathematics  Faculties  of  universities, 

for  he  was  completely  aware  of  the  increasingly  important  role 
played  Ly  present-day  mathematical  methods  in  solving  the 
problems  of  wave  diffraction.  As  an  active-proponent  of  more 
mathematical  training  for, -scientists  concerned  with  wave 
diffraction  theory,  Malyuzhinets  spared  no  efforts  to  master  for 
himself  new  mathematical  mc’hods  and  concepts.  ,  .  . 
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